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Abstract 

P-stars are compact stars made of up and down quarks in /3-equilibrium with 
electrons in a chromomagnetic condensate. P-stars are able to account for compact 
stars like RXJ 1856.5-3754 and RXJ 07204-3125, stars with radius comparable with 
canonical neutron stars, as well as super massive compact objects like SgrA*. We 
discuss p-stars endowed with super strong dipolar magnetic field which, following 
consolidated tradition in literature, are referred to as magnetars. We show that 
soft gamma-ray repeaters and anomalous X-ray pulsars can be understood within 
our theory. We find a well defined criterion to distinguish rotation powered pulsars 
from magnetic powered pulsars. We show that glitches, that in our magnetars 
are triggered by magnetic dissipative effects in the inner core, explain both the 
quiescent emission and bursts in soft gamma-ray repeaters and anomalous X-ray 
pulsars. We are able to account for the braking glitch from SGR 1900+14 and 
the normal glitch from AXP IE 2259+586 following a giant burst. We discuss and 
explain the observed anti correlation between hardness ratio and intensity. Within 
our magnetar theory we are able to account quantitatively for light curves for both 
gamma-ray repeaters and anomalous X-ray pulsars. In particular we explain the 
puzzling light curve after the June 18, 2002 giant burst from AXP IE 2259+586, so 
that we feel this last event as the Rosetta Stone for our theory. Finally, in Appendix 
we discuss the origin of the soft emission not only for soft gamma-ray repeaters and 
anomalous X-ray pulsars, but also for isolated X-ray pulsars. We also offer an 
explanation for the puzzling spectral features in IE 1207.4-5209. 
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1 INTRODUCTION 



In few years since their discovery pQ, pulsars have been identified with rotating neutron 
stars, first predicted theoretically by W. Baade and F. Zwicky [2|, endowed with a strong 
magnetic field |3J Ej. The exact mechanism by which a pulsar radiates the energy ob- 
served as radio pulses is still a subject of vigorous debate jSlEj, nevertheless the accepted 
standard model based on the picture of a rotating magnetic dipole has been developed 
since long time |H| • 

Nowadays, no one doubts that pulsars are neutron stars, even though it should be remem- 
bered that there may be other alternative explanations for pulsars. Up to present time it 
seems that there are no alternative models able to provide as satisfactory an explanation 
for the wide variety of pulsar phenomena as those built around the rotating neutron star 
model. However, quite recently we have proposed |3J a new class of compact stars, named 
p-stars, which is challenging the two pillars of modern astrophysics, namely neutron stars 
and black holes. We are, however, aware that such a dramatic change in the standard 
paradigm of relativistic astrophysics which is based on neutron stars and black holes needs 
a careful comparison with the huge amount of observations collected so far for pulsar and 
black hole candidates. In our opinion we feel that there are already clear observational 
evidences pointing towards the need of a drastic revision of the standard paradigm. So 
that, before addressing the main subject of the present paper, it is worthwhile to briefly 
discuss some observational evidences for p-stars and against neutron stars and black hole 
candidates. 

As concern black holes, we point out that the most interesting and intriguing aspect of 
our theory is that p-stars are able to overcome the gravitational collapse even for mass 
much greater 10 6 Mq. Indeed, from the equation of of state of degenerate up and down 
quarks in a chromomagnetic condensate described for the first time in Ref. we have 
on dimensional ground: 

where G is the gravitational constant, e c = e c /(gH) 2 , and e c is the central density. As a 
consequence we get: 

2 G M f(e c ) , , , i \ 

From Equation (jl.lj) we see that by decreasing the strength of the chromomagnetic con- 
densate we increase the mass and radius of the star. However, the ratio 2 G R M depends 
on e c only. It turns out that the function h(x) defined in Eq. (jl.2j) is less than 1 for 
any allowed values of e c [SJ. Thus, we infer that our p-stars do not admit the existence 
of an upper limit to the mass of a completely degenerate configuration. In other words, 
our peculiar equation of state of degenerate up and down quarks in a chromomagnetic 
condensate allows the existence of finite equilibrium states for stars of arbitrary mass. 
The accepted arguments for the evidence of black holes is based on the fact that there 
is spectroscopic evidence of compact objects with mass exceeding 3Mq. Indeed, a white 
dwarf cannot have a mass exceeding the Chandrasekhar limit, about 1.4 Mq, while even 
for neutron stars there is a maximum mass which probably is about 3 Mq . Then, com- 
pact objects with mass exceeding 3Mq are classified as black holes. However, as we 
argued before, this argument cannot distinguish between a massive p-star and a black 
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hole. Indeed, the fundamental difference between massive p-stars and black holes resides 
in the lack of stellar surface in black holes. As it is well known, from general relativity 
it follows that the black hole boundary is a geometric surface called the event horizon. 
Recently, there are claims in literature for evidence of event horizons. For instance, in 
Ref. [12] it is claimed that the X-ray luminosities of black hole candidates in quiescence 
are much less than the corresponding X-ray luminosities of compact solar mass stellar 
systems. However, the authors of Ref. [13] argued that it is not correct to compare only 
the X-ray luminosities. If one compares the bolometric luminosities, then it turns out 
that there are no observable differences between compact solar mass stellar systems and 
black holes candidates. So that, up to now there is no compelling evidence in favour 
of event horizons. On the other hand, interestingly enough the authors of Ref. [H], by 
using the standard analysis of magnetic propeller effect jTH] for pulsar in low mass X- 
ray binaries, found that the spectral properties of galactic black hole candidates could 
be accounted for by compact objects with an intrinsic magnetic moment. Subsequently, 
in Ref. [TH] these authors, extending their analysis to active galactic nuclei, showed how 
a standard accretion disk can interact with the intrinsically magnetized central compact 
object to drive low state jets. Even though these authors believe that massive intrinsically 
magnetized central objects can be accounted for within general relativity as highly red 
shifted, extremely long lived, collapsing, radiating objects [T7|[THj, it is evident that mas- 
sive p-stars endowed with magnetic fields are indeed the natural candidates for massive 
compact objects with an intrinsic magnetic moment. 

As we will discuss in a future paper, in p-stars there is a natural mechanism to generate a 
dipolar magnetic field. As a matter of fact, it turns out that the generation of the dipolar 
magnetic field is enforced by the formation of a dense inner core composed mainly by 
down quarks. As a consequence the surface dipolar magnetic field Bs is proportional to 
the strength of the chromomagnetic condensate gH. More precisely we have (here and in 
the following we shall adopt natural units fi = c = ks = 1): 



where e is the electric charge, R and R c are the stellar and inner core radii respectively. 
It is interesting to observe that for p-stars with canonical mass M ~ 1.4 Mq we get 



y/gH ~ 0.55 GeV. On the other hand, massive p-stars with M ~ 10 Mq require a 

chromomagnetic condensate gH smaller by about a factor 4 — 5 with respect to canonical 
p-stars. So that, from Eq. (|1.3|) it follows that massive p-stars are characterized by sur- 
face magnetic fields reduced by less than one order of magnitude with respect to pulsar 
magnetic fields. 

In addition, as we said, there are finite equilibrium states for stars of arbitrary mass. For 
instance, SgrA*, the super massive compact object at the galactic center [TH], could be 
interpreted as a p-star with mass M ~ 3.2 10 6 Mq and radius R ~ 1.4 10 7 Km. The cor- 
responding strength of the chromomagnetic condensate turns out to be \fg~H ~ 0.4 MeV. 
Recent CHANDRA observations [2*0*] |2*T] of diffuse emission around the galactic center 
have confirmed that SgrA* is underluminous in X-ray by a factor of about 10 7 — 10 8 
compared to the standard thin accretion disk model. The low luminosity of SgrA* may 
be explained within the standard paradigm either by accretion at a rate far below the 
estimate Bondi rate, or accretion at the Bondi rate of gas that is radiating very inef- 
ficiently. Since the gas from the winds of the surrounding young massive stars should 
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be able to maintain the accretion rate at a sizeable fraction of the Bondi accretion rate, 
which has been estimate [20] to be about 10~ 6 Mq years' 1 , it is believed that only the 
latter possibility is viable. However, in Ref. [22], using close stars as probe of accretion 
flow in SgrA*, it has been pointed out that non radiative accretion flows are constrained 
to accretion rates no larger than 10~ 7 Mq year s~ l . It is not yet clear if this constrain 
could be reconciled with the accretion rate estimate in Ref. [20] • Furthermore, Zhao et 
al. [23] reported the presence of 106 days cycle variability at centimeter wavelengths in 
the radio flux density of SgrA* . This peculiar periodicity has been confirmed by observa- 
tions at millimeter wavelengths [21]. The very low X-ray luminosity and the periodicity 
in the flux density of SgrA* look puzzling within the standard interpretation based on 
accreting black holes, while these can be accounted for if we assume that SgrA* is a super 
massive p-star. Indeed, the periodicity is naturally explained assuming a rotation period 
P ~ 106 days ~ 9.2 10 6 sec. Moreover, from the strength of the chromomagnetic con- 
densate and from Eq. (jl.Hj) we estimate that the dipolar surface magnetic field of SgrA* 
is reduced by about 10~ 6 with respect to pulsar magnetic fields. So that Bs gr A* should 
lie in the range 10 6 — 10 9 Gauss. From Bs gr A* % 10 9 Gauss, we infer for the age of 
SgrA* t ~ 10 10 years, i.e. SgrA* is a primordial p-star. Finally, the low X-ray quiescent 
luminosity could be interpreted as thermal emission from the stellar surface. 
As discussed in Refs. ED] p-stars are compact stars made of up and down quarks in 
/3-equilibrium with electrons in an abelian chromomagnetic condensate. It turns out that 
these compact stars are more stable than both neutron stars and strange stars whatever 
the value of the chromomagnetic condensate y 'gH. In other words, p-stars, once formed, 
are absolutely stable. The logical consequence is that now we must admit that supernova 
explosions give rise to p-stars. In other words, we are lead to identify pulsars with p-stars 
instead of neutron or strange stars. Such a dramatic change in the standard paradigm 
of relativistic astrophysics has been already advanced in our previous paper [§| where we 
suggested that, if we assume that pulsars are p-stars, then we could solve the supernova 
explosion problem. As is well known, the binding energy is the energy released when the 
core of an evolved massive star collapses. Actually, only about one percent of the energy 
appears as kinetic energy in the supernova explosion [TT|. Now, in Ref. [§] we showed that 
there is an extra gain in kinetic energy of about 1 — 10 foe {I foe = 10 51 erg), which 
is enough to fire the supernova explosions. Further support to our theory comes from 
cooling properties of p-stars. In fact, we found that cooling curves of p-stars compare 
rather well with available observational data. 

In our previous papers EI] we showed that p-stars are also able to account for compact 
stars with R < 6 Km. In particular, we convincingly argued that the nearest isolated pul- 
sars RXJ 1856.5-3754 and RXJ 0720.4-3125 (for a recent review see j2SH2E]) are p-stars. 
From the X-ray emission spectrum we argued that the most realistic interpretation is that 
these objects are compact p-stars with M ~ 0.8 Mq and R ~ 5 Km. However, it should 
be stressed that in the observed spectrum there is also an optical emission in excess over 
the extrapolated X-ray blackbody. By interpreting the optical emission as a Rayleigh- 
Jeans tail of a thermal blackbody emission, one finds that the optical data are also fitted 
by the blackbody model yielding an effective radius R°° > 16 Km 12 Q pc [2*7] . However, 
interestingly enough, quite recently the distance measurement of RXJ 1856.5-3754 has 
been reassessed and it is now estimated to be at 180 pc instead of 120 pc [21]. This new 
determination of the distance of RXJ 1856.5-3754 rules out the two blackbody interpre- 
tation of the spectrum, for this model leads now to an effective radius R°° > 24 Km, 



3 



which is too large for a neutron star. Thus, the new determination of the distance of 
RXJ 1856.5-3754 strongly supports our p-star theory, and indicates clearly that the faint 
optical emission originates in the magnetosphere. In general, it remarkable that isolated 
X-ray pulsars do display a faint soft emission, in excess over the extrapolated X-ray 
thermal emission, which is best fitted by a non thermal power law. The origin of this 
faint emission is puzzling. Nevertheless, our previous considerations point toward a gen- 
eral mechanism in the magnetosphere responsible for the faint emission. This problem, 
which to the best of our knowledge has never addressed before, is thoroughly analyzed 
in Appendix, where we show that a subtle quantum mechanical effect related to strong 
enough magnetic fields in the polar cap regions leads to a faint non thermal power law 
soft emission. 

Quite recently it has been proposed that the compact accreting object in the famous X-ray 
binary Herculses X-l is a strange star [22] • This proposal was based on the compari- 
son of a phenomenological mass-radius relation for Herculses X-l jHHj with theoretical 
M — R curves for neutron and strange stars. The analysis in Ref. [2H] has, however, 
been criticized by the authors of Ref. j^J. These authors, using a new mass estimate 
together with a revised distance, which leads to a somewhat higher X-ray luminosity, 
argued that the hypothesis that Herculses X-l is a neutron star is not disproved. As 
a matter of fact, the authors of Ref. [31] found that there is marginal consistency with 
observations if one adopts for neutron stars a very soft equation of state. At the same 
time, these authors pointed out that the hypothesis of a strange star can be ruled out 
since the theoretical curves no longer intercept the observational relations within the per- 
mitted mass range. Recent observations of millisecond pulsars in the globular cluster 
Terzan 5 using the Green Bank Telescope |32] indicated that al least one of the pulsar is 
more massive than M ~ 1.7 Mq. Even more, there is emerging observational evidences 
for pulsars with mass in excess of 1.6 Mq. For instance, the pulsar in Vela X-l has mass 
M = 1.86 ± 0.16 Mq [53] I34j . The very existence of such massive pulsars constrains the 
equation of state of matter in neutron stars. In fact it seems that soft equations of state 
are almost certainly ruled out. As a consequence we infer that the compact accreting 
pulsar in Herculses X-l cannot be a strange star nor a neutron star. On the other hand, 
theoretical M — R curves for p-stars are compatible with the phenomenological mass- 
radius relation for Herculses X-l. Indeed, we find that the pulsar in Herculses X-l 
could be described by a p-star with M ~ 1.5 Mq, R ~ 6.5 Km, and y 'gH ~ 0.62 GeV. 
In the present paper we investigate the properties of p-stars with super strong surface 
magnetic field. As we shall show, these p-stars allow us reach a complete understand- 
ing of several puzzling observational aspects of anomalous X-ray pulsars (AXPs) and 
soft gamma-ray repeaters (SGRs). Anomalous X-ray pulsars and soft gamma-ray re- 
peaters are two class of intriguing objects that in our opinion are challenging the stan- 
dard paradigm based on neutron stars. For a recent review on the observational proper- 
ties of anomalous X-ray pulsars see Refs. |3Sl EE1 EZ[ , for soft gamma-ray repeaters see 
Refs. jSHl EH] • Recently, these two groups have been linked by the discovery of persistent 
emission from soft gamma ray repeaters that is very similar to anomalous X-ray pulsars, 
and bursting activity in anomalous X-ray pulsars quite similar to soft gamma ray re- 
peaters (see, for instance Refs. [471] l4*T]). 

Duncan and Thompson [12] and Paczyhski [13] have proposed that soft gamma-ray re- 
peaters are pulsars whose surface magnetic fields exceed the QED critical magnetic field: 
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Indeed, Duncan and Thompson in Ref. [12] refer to these pulsar as magnetars. In partic- 
ular Duncan and Thompson 03] argued that the soft gamma-ray repeater bursts and 
quiescent emission were powered by the decay of an ultra-high magnetic field. This inter- 
pretation is based on the observations that showed that these peculiar pulsars are slowing 
down rapidly, with an inferred magnetic dipole field much greater than the quantum crit- 
ical field Bqed-i while producing steady emission at a rate far in excess of the rotational 
kinetic energy loss. The identification of anomalous X-ray pulsars with magnetars was 
more recently supported by the similarity of anomalous X-ray pulsar emission to that of 
soft gamma-ray repeaters in quiescence. This was confirmed by the detection of SGR-like 
bursts from two anomalous X-ray pulsars jJOJ . 

In the standard neutron star theory, magnetars ought to be born with millisecond initial 
period to ensure vigorous dynamo process to occur This mechanism should generate 
huge surface dipolar magnetic field up to 10 15 Gauss, and even stronger interior fields. 
However, strong magnetic fields in excess of the critical field Bqed would squeeze electrons 
into the lowest Landau levels. In this conditions, the electron gas pressure transverse to 
the magnetic field may vanish leading to a transverse collapse of the star jI7j. On the other 
hand, p-stars do not share this stability problem. In fact, Eq. ()1.3|) shows that the dipolar 
magnetic field in p-star is a tiny effect with respect to the chromomagnetic condensate. 
Moreover, the stability of p-stars is due to the quark pressure, while the electron pressure 
is almost completely negligible. Indeed, from Eq. (J1.3j) it follows that canonical p-stars 
could support dipolar magnetic fields up to 10 16 Gauss. As we said before, in standard 
magnetars appropriate conditions for true dynamo mechanism could exist if neutron star 
is born with a very short period. However, up to now there is no direct observational 
evidences for such short initial periods in radio pulsars or in young supernova remnants. 
On the contrary, the peculiar mechanism to generate dipolar magnetic fields in p-stars 
indicates that huge magnetic fields require rather large initial period. It is worthwhile to 
stress that our mechanism for the generation of dipolar magnetic fields in pulsars solves 
in a natural way the puzzling discrepancy between characteristic ages and true ages which 
is displayed by at least two anomalous X-ray pulsars. 
As it is well known, the pulsar characteristic ages is defined as: 



while the true age is given by: 



The true age can be significantly smaller than r c if the initial period Pq is close to P. 
Both anomalous X-ray pulsars IE 1841-045 and IE 2259+586 have r c greater than the 
ages estimated from their respective supernova remnants (see Table 1 in Ref. 36J). In 
particular, for IE 2259+586 the characteristic age r c is more than one order of magnitude 
larger than the age of SNR G109. 1-0.1. We may solve this large discrepancy if we assume 
the initial period Pq ~ 6.93 sec. In the case of IE 1841-045 we find that P ~ 10.6 sec 




(1.6) 
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reconciles the age discrepancy. Note that this solution cannot be adopted within the stan- 
dard neutron star theory, for huge magnetic fields can be generated only for stars born 
with very small initial period. On the other hand, the peculiar mechanism to generate 
dipolar magnetic fields in p-stars favors initial strong magnetic fields in slowly rotating 
stars. In the next Section we shall present further phenomenological evidences in support 
of our point of view. 

The main purpose of this paper is to discuss in details p-stars endowed with super strong 
dipolar magnetic field which, following well consolidated tradition in literature, will be re- 
ferred to as magnetars. In particular we will show that, indeed, soft gamma-ray repeaters 
and anomalous X-ray pulsars can be understood within our theory. Whenever possibly, 
we shall critically compare our theory with the standard paradigm based on neutron stars. 
The plan of the paper is as follows. In SectionEJwe discuss the phenomenological evidence 
for the dependence of pulsar magnetic fields on the rotational period. We argue that there 
is a well defined criterion which allows us to distingue between rotation powered pulsars 
and magnetic powered pulsars. We explicitly explain why the recently discovered high 
magnetic field radio pulsars are indeed rotation powered pulsars. In Section |3] we intro- 
duce the radio death line, which in the P — P plane separated radio pulsars from radio 
quiet magnetic powered pulsars, and compare with available observational data. Section 0] 
is devoted to the glitch mechanism in our magnetars and their observational signatures. 
In Section |4~T1 we compare glitches in SGR 1900+14 and IE 2259+586, our prototypes for 
soft gamma-ray repeater and anomalous X-ray pulsar respectively. Sections 14.21 and 14.31 
are devoted to explain the origin of the quiescent luminosity, the bursts activity and the 
emission spectrum during bursts. In Section l4~4l we discuss the puzzling anti correlation 
between hardness ratio and intensity. In Section El we develop a general formalism to cope 
with light curves for both giant and intermediate bursts. In Sections 15.11 through 15.41 we 
careful compare our theory with the available light curves in literature. In particular, we 
are able to account for the peculiar light curve following the June 18, 2002 giant burst 
from the anomalous X-ray pulsar IE 2259+586. Finally, we draw our conclusions in Sec- 
tion El In Appendix we face with the problem of the origin of soft emission in isolated 
X-ray pulsars. We argue that the soft emission originates in the magnetosphere and it 
can be ascribed to a subtle quantum mechanical effect related to strong magnetic fields 
in the polar cap regions. 

2 ROTATION VERSUS MAGNETIC POWERED PULSARS 

In Ref. jH] we introduced p-stars, a new class of compact quark stars made of almost mass- 
less deconfined up and down quarks immersed in a chromomagnetic field in /3-equilibrium. 
The structure of p-stars is determined once the equation of state appropriate for the de- 
scription of deconfined quarks and gluons in a chromomagnetic condensate is specified. 
In particular, the chemical potentials satisfy the constrains arising from f3- equilibrium 
and charge neutrality: 
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Figure 1: Inferred magnetic field B$ plotted versus stellar period for 1194 pulsars taken 
from the ATNF Pulsar Catalog jlHl- Full red line corresponds to Eq. ()2.5|) with B\ ~ 
1.3 10 13 Gauss. 



where: 
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gHfx u , n d 
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From previous equations it follows that fid > fi u > fi e - Moreover, it turns out that 
the chemical potentials are monotonic decreasing smooth functions of the distance from 
the center of the star. In general, the quark chemical potentials fid and fi u are smaller 
that the strength of the chromomagnetic condensate y/gH. So that, up and down quarks 
occupy the lowest Landau levels. However, for certain values of the central energy density 
it happens that -j= > 1 in the stellar core. Thus, a fraction of down quarks must jump 



into higher Landau levels leading to a central core with energy density e c somewhat greater 
than the energy density outside the core. Now, these quarks in the inner core produce 
a vector current in response to the chromomagnetic condensate. Obviously, the quark 
current tends to screen the chromomagnetic condensate by a very tiny amount. However, 



since the down quark has an electric charge 



e, the quark current generates in the 



core a uniform magnetic field parallel to the chromomagnetic condensate with strength: 



B c ~ 
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(2.4) 



Outside the core the magnetic field is dipolar leading to surface magnetic field given 
by Eq. p. 3)1 . In general the formation of the inner core denser than the outer core is 
contrasted by the centrifugal force produced by stellar rotation. Since the centrifugal 
force is proportional to the square of the stellar rotation frequency, this leads us to argue 
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that the surface magnetic field strength is proportional to the square of the stellar period: 



Bs - B x (JL) , (2.5) 
\1 sec J 

where B\ is the surface magnetic field for pulsars with nominal period P = 1 sec. Indeed, 
in Fig. ^ we we display the surface magnetic field strength Bs inferred from (for instance, 
see Ref. @E]): 

B s ~ 3.1 10 19 v r PP Gauss , (2.6) 

versus the period. Remarkably, assuming Bi ~ 1.3 10 13 Gauss, we find the Eq. (12. 5|) ac- 
counts rather well the inferred magnetic field for pulsars ranging from millisecond pulsars 
up to anomalous X-ray pulsars and soft-gamma repeaters. As a consequence of Eq. ()2.5j) . 
we see that the dipolar magnetic field is time dependent. In fact, it is easy to find: 



B s (t) ~ B 1 + 2 - t) , (2.7) 




where B indicates the magnetic field at the initial observation time. Note that Eq. ()2.7|) 
implies that the magnetic field varies on a time scale given by the characteristic age. Equa- 
tion ()2.7|) leads to remarkable consequences discussed in Ref. jSHl- Indeed, in Ref. |5U] . 
starting from Equation ([2.7)1 . we discussed a general mechanism which allows to explain 
naturally both radio and high energy emission from pulsars. We also discuss the plasma 
distribution in the region surrounding the pulsar, the pulsar wind and the formation of 
jet along the magnetic axis. We also suggested a plausible mechanism to generate pulsar 
proper motion velocities. In particular, in our emission mechanism there is a well de- 
fined geometric mapping between frequency and distance from the star which seems to be 
consistent with observations. In particular, in the recently detected binary radio pulsar 
system J0737-3039 jSH E2] it has been reported jHHJ the detection of features similar to 
drifting subpulses with a fluctuation frequency which is exactly the beat frequency be- 
tween the period of the two pulsars. This direct influence of the electromagnetic radiation 
from one pulsar on the electromagnetic emission from the other pulsar can be naturally 
accounted for within our emission mechanism due to the geometric mapping between 
emission frequencies and distances. On the other hand, that effect cannot easily recon- 
ciled with the generally accepted model for pulsar radio emission which involves coherent 
radiation from very energetic e + — er pairs. 

It is widely accepted that pulsar radio emission is powered by the rotational energy: 

Er = \ l uj 2 , (2.8) 
so that, the spin-down power output is given by: 

— Er — — I lo Co = 4 7T 2 7 . (2.9) 

On the other hand, an important source of energy is provided by the magnetic field. 
Indeed, the classical energy stored into the magnetic field is: 



E B = ^- / d 3 rB 2 (r) , (2.10) 
8vr J r>R 
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Assuming a dipolar magnetic field: 
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B(r) = B s I — ) for r > R , (2.11) 



Eq. (I2~TU1) leads to: 

E B = I B 2 S R 3 . (2.12) 
o 

Now, from Eq. (|2.7jl the surface magnetic field is time dependent. So that, the magnetic 
power output is given by: 

E b = I Bl R 3 *p . (2.13) 

For rotation-powered pulsars it turns out that E B <C — E R . However, if the dipolar 
magnetic field is strong enough, then the magnetic power Eq. (|2.13|) can be of the order, or 
even greater than the spin-down power. Thus, we may formulate a well defined criterion 
to distinguish rotation powered pulsars from magnetic powered pulsars. Indeed, until 
Eb < — Ejt there is enough rotation power to sustain the pulsar emission. On the other 
hand, when Eg > — Er all the rotation energy is stored into the increasing magnetic 
field and the pulsar emission is turned off. In fact, in the next Section we will derive the 
radio death line, which is the line that in the P — P plane separates rotation-powered 
pulsars from magnetic-powered pulsars. In the remainder of this Section, we discuss the 
recently detected radio pulsar with very strong surface magnetic field. We focus on the 
two radio pulsars with the strongest surface magnetic field: PSR J1718-3718 and PSR 
J1847-0130 . These pulsar have inferred surface magnetic fields well above the quantum 
critical field Bq ED above which some models [51] predict that radio emission should not 
occur. In particular, we have: 

PSRJ1718 - 3718 E2 P ~ 3.4 sec , Be ~ 7.4 10 13 Gauss 



s . . _ __ , 
PSRJ1847 - 0130 [SHI P - 6.7 sec , B s ~ 9.4 10 13 Gauss . K ' J 

Both pulsars have average radio luminosities and surface magnetic fields larger than that 
of AXP IE 2259+586. Now, using Eqs. flTOD , Q , together with I =\MR 2 , and 
Eq. (I2~T1) we get: 

PSRJ1718 - 3718 - E R ~ 3.4 10 45 erg — , Eg ~ 4.7 10 44 erg - , 

F . F . (2.15) 

P P 
PSRJ18A7 - 0130 - E R ~ 8.8 10 44 erg — , E B ~ 7.5 10 44 erg — . 

We see that in any case: E B < —E R , so that there is enough rotational energy to power 
the pulsar emission. 



3 RADIO DEATH LINE 

As discussed in previous Section, until E B < — Er the rotation power loss sustains the 
pulsar emission. We have already shown that this explain the otherwise puzzling pulsar 
activity for pulsars with inferred magnetic fields well above the critical field Bq B d- Even 
more, the two pulsars PSR J 17 18-37 18 and PSR J 1847-0130 have magnetic fields which 
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are well above the magnetic field of the anomalous X-ray pulsar AXP IE 2259+586. Up 
to now, it was unclear how these high-field pulsars and anomalous X-ray pulsars can 
have such similar spin-down parameters but vastly different emission properties. We have 
offered a natural explanation for the pulsar activity for high-field radio pulsars. In this 
Section we explain why anomalous X-ray pulsars and soft gamma repeaters are radio 
quiet pulsars. When Eb > — Er all the rotation energy is stored into the increasing 
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Figure 2: Period derivative plotted versus stellar period for 1194 pulsars taken from the 
ATNF Pulsar Catalog fH]. Full blue line corresponds to Eq. (jTOj) 

magnetic field and the pulsar emission is turned off. As a consequence pulsars with strong 
enough magnetic fields are radio quiet. Accordingly we see that the condition: 

E B = -E R . (3.1) 

is able to distinguish rotation powered pulsars from magnetic powered pulsars. Now, 
using |48 j 

we recast Eq. ()3.1j) into: 

P 3 P = 16tt 4 P 3 . (3.3) 
Using the canonical radius R ~ 10 Km, we get: 

3 log(P) + log(P) ~ - 10 . (3.4) 

Equation ()3.4|) is a straight line, plotted in Fig.|2l in the log(P)— log(P) plane. In Figure. 121 
we have also displayed 1194 pulsars taken from the ATNF Pulsar Catalog |49] . We see 
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that rotation powered pulsars, ranging from millisecond pulsars up to radio pulsars, do 
indeed lie below our Eq. (|3.4|) . Note that in Fig. |2]the recently detected high magnetic 
field pulsars are not included. However, we have already argued in previous Section that 
these pulsars have spin parameters which indicate that these pulsars are rotation powered. 
On the other hand, Fig. |21 shows that all soft gamma-ray repeaters and anomalous X-ray 
pulsars in the ATNF Pulsar Catalog lie above our line Eq. ()3.4j) . In particular, in Fig. |21 
the pulsar above and nearest to the line Eq. ()3.4j) corresponds to AXP IE 2259+586. So 
that, we see that our radio dead line, Eq. ()3.4|) . correctly predicts that AXP IE 2259+586 
is not a radio pulsar even though the magnetic field is lower than that in radio pulsars 
PSR J1718-3718 and PSR J1847-0130. We may conclude that pulsars above our dead 
line are magnetars, i.e. magnetic powered pulsars. The emission properties of magnetars 
are quite different from rotation powered pulsars. The emission from magnetars consists 
in thermal blackbody radiation form the surface. In addition, it could eventually also be 
present a faint power-law emission superimposed to the thermal radiation. As discussed 
in Appendix, this soft faint emission is caused by the thermal radiation reprocessed in the 
magnetosphere. In Ref. [TU] we suggested that RXJ 1856.5-3754 is exactly in this state. 
On the other hand, the energy stored into the magnetic field can be released if the star 
undergoes a glitch. Indeed, as thoroughly discussed in the next Section, glitches originate 
from dissipative effects in the inner core of the star leading to a decrease of the strength 
of the core magnetic field. So that, soon after the glitch there is a release of magnetic 
energy. We have already suggested in Ref. [TU] that this picture is consistent with the 
long-term variability in the X-ray emission of RXJ 0720.4-3125. Remarkably, a recent 
timing analysis of the isolated pulsar RXJ 0720.4-3125 performed in Ref. |HZ| suggested 
that, among different possibilities, glitching may have occurred in this pulsar. 

4 GLITCHES IN MAGNETARS 

In p-stars there is a natural mechanism to generate a dipolar magnetic field, namely the 
generation of the dipolar magnetic field is enforced by the formation of a dense inner core 
composed mainly by down quarks. A quite straightforward calculation, which will be 
presented elsewhere, leads to the conclusion that down quarks in the inner core produce a 
vector current in response to the chromomagnetic condensate. This quark current, in turn, 
generates in the core a uniform magnetic field parallel to the chromomagnetic condensate 
with strength given by Eq. (|2.4j) . Outside the core the magnetic field is dipolar. Now, we 
note that the inner core is characterized by huge conductivity, while outer core quarks 
are freezed into the lowest Landau levels. So that, due to the energy gap between the 
lowest Landau levels and the higher ones, the quarks outside the core cannot support any 
electrical current. As a consequence, the magnetic field in the region outside the core is 
not screened leading to our previous Eq. (jl.3|) . For later convenience, after taking into 
account Eq. (|2.4j) . we rewrite Eq. (jl.3j) as: 

B S * B< ,B^^- 9 H. (4.1, 

In Figure El we display a schematic view of the interior of a p-star. 

The presence of the inner core uniformly magnetized leads to well defined glitch mech- 
anism in p-stars. Indeed, dissipative effects, which are more pronounced in young stars, 
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Figure 3: Schematic view of the interior of a p-star. R c and R represent the inner core 
and stellar radii respectively; e c is the energy density of the inner core, e is the energy 
density outside the core. 



tend to decrease the strength of the core magnetic field. On the other hand, when B c de- 
creases due to dissipation effects, then the magnetic flux locally decreases and, according 
to Lenz's law, induces a current which resists to the reduction of the magnetic flux. This 
means that some quarks must flow into the core by jumping onto higher Landau levels. 
In other words, the core radius must increase. Moreover, due to very high conductivity 
of quarks in the core, we have: 

B c R\ ~ constant , (4.2) 

which implies: 

SB ' -2^*0, (4.3) 



B c R c 

^ ^ - i ^ . (4.4) 
R c 2 B c 

Equation (J4.4j) confirms that to the decrease of the core magnetic field, 5B C < 0, it 
corresponds an increase of the inner core radius 5R C > 0. This sudden variation of the 
radius of the inner core leads to glitches. Indeed, it is straightforward to show that the 
magnetic moment: 

m = B s R 3 = B c Rl (4.5) 
where we used Eq. (|4.1|) . must increase in the glitch. Using Eq. (|4.4jl . we get: 

5 m 8B C 8R C 5R C 

= ~5~ + 3 ~^~~^ >0 ' ( 4 - 6 ) 
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Another interesting consequence of the glitch is that the stellar radius R must decrease, 
i.e. the star contracts. This is an inevitable consequence of the increase of the inner core, 
which is characterized by an energy density e c higher then the outer core density e. As 
a consequence the variation of radius is negative: 5R < (see Fig. El). In radio pulsar, 
where the magnetic energy can be neglected, conservation of the mass leads to: 



5R £ c — £ ( R c \ 3 8 R c 



R £ \ R I R, 



(4.7) 



where e is the average energy density. In general, we may assume that ^=-^ is a constant 
of order unity. So that Eq. ()4.7|) becomes: 



5R {R r \ 3 5R r 



R \ R J R r 

,3 



(4.8) 



Note that the ratio (jf) can be estimate from Eq. (j4.1j) once the surface magnetic field 

is known. We find that, even for magnetars, (-7^) 3 is of order 10~ 2 or less. So that our 
Eqs. (@~U) and (@~BJ) show that: 

sr < , - 5 4« 5 -^- < 49 > 

it it c 

As is well known, because the external magnetic braking torque, pulsars slow down ac- 
cording to (e.g. see Ref. j4*%|): 

- oc m 2 r 1 z/ 3 . (4.10) 

So that: 

^ = 2^-^ + 3^. (4.11) 
v ml v 



From conservation of angular momentum we have: 

lyl - 1—1 • 

Moreover, from observational data it turns out that: 



(4.12) 



, 80. ,5u. , 
l-l » l-l , (4-13) 



so that Eq. ()4.11j) becomes: 



50 n 5m n 5 R c 5 B c ,,+,s 

— ~ 2 — ~ 2 ~ , (4.14) 

m R c B c 

where we used Eqs. ()4.6|) and ()4.4|) . Equation f)4.14j) does show that the variation of the 
radius of the inner core leads to a glitch. 

In rotation powered pulsar, starting from Eq. (|4.7|) one can show that |^| ~ 1^1- So that, 
taking into account Eq. (|4.9jl we recover the phenomenological relation Eq. (|4.13)1 . A full 
account of glitches in radio pulsar will be presented elsewhere. Glitches in magnetars are 
considered in the next Section, where we show that, indeed, Eqs. ()4.9|) and (|4.13J) hold 
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even in magnet ars. 

The most dramatic effect induced by glitches in magnetars is the release of a huge amount 
of magnetic energy in the interior of the star and into the magnetosphere. To see this, let 
us consider the energy stored into the magnetic field in the interior of the magnetar. We 
have: 



rpint 



6 Rl Bl 



1 



i7T 



R 



— / d 6 r 



R c 



D, [ ^ 

r 



-I 2 



(4.15) 



where the first term is the energy stored into the core where the magnetic field is uniform. 
The variation of the magnetic energy Eq. (|4.15jl caused by a glitch is easily evaluated. 
Taking into account Eq. (|4.4|) and (-7^) 3 <C 1, we get: 



6 Eg* 



SRc _ 1 ^3 B 2 5R c 

R r 2 R r 



1 R 3 
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SRc 

R r 



(4.16) 

Equation ()4.16|) shows that there is a decrease of the magnetic energy. So that after a 
glitch in magnetars a huge magnetic energy is released in the interior of the star. We 
shall see that this energy is enough to sustain the quiescent emission. On the other 
hand, the glitch induces also a sudden variation of the magnetic energy stored into the 
magnetosphere. Indeed, from Eq. (J4.1j) we find: 



SB, 
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(4.17) 



Thus, the magnetic energy stored into the magnetosphere: 
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(4.18) 



increases by: 



5E e B xt 



- R 3 B\ 

3 5 



B s 2 R J 3 



SJh 
Bs 



> 



(4.19) 



This magnetic energy is directly injected into the magnetosphere, where it is dissipated by 
well defined physical mechanism discussed in Section 14. 3^ and it is responsible for bursts 
in soft gamma-ray repeaters and anomalous X-ray pulsars. 

To summarize, in this Section we have found that dissipative phenomena in the inner core 
of a p-star tend to decrease the strength of the core magnetic field. This, in turn, results 
in an increase of the radius of the core SR C > 0, and in a contraction of the surface of the 
star, 5R < 0. We have also show that the glitch releases an amount of magnetic energy 
in the interior of the star and injects magnetic energy into the magnetosphere, where it 
is completely dissipated. Below we will show that these magnetic glitches are responsible 
for the quiescent emission and bursts in gamma-ray repeaters and anomalous X-ray pul- 
sars. Interestingly enough, in Ref. [SH] it was shown that SGR events and earthquakes 
share several distinctive statistical properties, namely: power-law energy distributions, 
log-symmetric waiting time distributions, strong correlations between waiting times of 
successive events, and weak correlations between waiting times and intensities. These 
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statistical properties of bursts can be easily understood if bursts originate by the release 
of a small amount of energy from a reservoir of stored energy. As a matter of fact, in 
our theory the burst activity is accounted for by the release of a tiny fraction of magnetic 
energy stored in magnetars. Even for giant bursts we find that the released energy is 
a few per cent of the magnetic energy. Moreover, the authors of Ref. [2H] noticed that 
there is a significantly statistical similarity between the bursts from SGR 1806-20 and the 
microglitches observed from the Vela pulsar with |^| ~ 1CT 9 . So that we see that these 
early statistical studies of bursts are in complete agreement with our theory for bursts in 
magnetars. Even more, we shall show that after a giant glitch there is an intense burst 
activity quite similar to the settling earthquakes following a strong earthquake. 

4.1 BRAKING GLITCHES 

In Section HI we found that magnetic glitches in p-stars lead to: 

60 SB ° >0. (4.20) 



v B c 

Since there is variation of both the inner core and the stellar radius, the moment of inertia 
of the star undergoes a variation of 51. It is easy to see that the increase of the inner core 
leads to an increase of the moment of inertia /; on the other hand, the reduction of the 
stellar radius implies 51 < 0. In radio pulsar, where, by neglecting the variation of the 
magnetic energy, the conservation of the mass leads to Eq. (J4.8)) . one can show that: 

51 5R . 

T " If < • (421) 



Moreover, from conservation of angular momentum: 



5J_ 5v_ 
T ~ V 



(4.22) 



it follows: 

5v 5R 5B C 5v 

< ~ = "If « " -t H T ■ (423) 

For magnetars, namely p-stars with super strong magnetic field, the variation of magnetic 
energy cannot be longer neglected. In this case, since the magnetic energy decreases, we 
have that the surface contraction in magnetars is smaller than in radio pulsars. That 
means that Eq. ()4.9j) holds even for magnetars. Moreover, since in radio pulsars we 
known that: 

— = ~~f - "~W~ • <««) 

V 1 it 

we see that in magnetars the following bound must hold: 



5_R 
~R 



< lO" 6 . (4.25) 



As a consequence we may write: 



51 {5I\ f5I\ f5I\ 5R 



/ \ I ) — „ V / / sur f \ I J sur f R 
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As we show in a moment, the variation of the moment of inertia induced by the core is 
positive. So that if the core contribution overcomes the surface contribution to SI we 
have a braking glitch where — ^ = ^ > 0. 

We believe that the most compelling evidence in support to our proposal comes from 
the anomalous X-ray pulsar AXP IE 2259+586. As reported in Ref. 60j, the timing 
data showed that a large glitch occurred in AXP IE 2259+586 coincident with the 2002 
June giant burst. Remarkably, at the time of the giant flare on 1998 August 27, the soft 
gamma ray repeater SGR 1900+14 displayed a discontinuous spin-down consistent with 
a braking glitch [61J. Our theory is able to explain why AXP IE 2259+586 displayed a 
normal glitch, while SGR 1900+14 suffered a braking glitch. To see this, we recall the 
spin-down parameters of these pulsars: 



SGR 1900+14 P ~ 5.16 sec, P ~ 1.1 10~ 10 , B s ~ 7.4 10 14 Gauss , 
AXP IE 2259+586 P ~ 6.98 sec , P ~ 2.0 10~ 14 , B s ~ 1.2 10 13 Gawss . 



(4.27) 



For canonical magnetars with M ~ 1.4 Mq and radius R ~ 10 Km, we have ^/gif ~ 
0.55 GeV. So that, using 1 GeV 2 ~ 5.12 10 19 Gawss, we rewrite Eq. fETT]) as: 

fi 5 ~ 1.54 10 16 (^] Gauss . (4.28) 



R 



Combining Eqs. (J4~2"T|) and (|4~2"%j) we get: 



SGR 1900+14 ( ^ ] ~ 4.81 10~ 2 
.4 AT ././:' 1259+586 ( ) - 0.7s 



(4.29) 



According to Eqs. (|4.20|) . ()4.22j) and (|4.26|1 . to evaluate the sudden variation of the fre- 
quency and frequency derivative, we need SR and 5R C . These parameters can be estimate 
from the energy released during the giant bursts. In the case of AXP IE 2259+586, the 
giant 2002 June burst followed an intense burst activity which lasted for about one year. 
The authors of Ref. [Hill, assuming a distance of 3 kpc to IE 2259+586, measured an 
energy release of 2.7 10 39 ergs and 2.1 10 41 ergs for the fast and slow decay intervals, re- 
spectively. Due to this uncertainty, we conservatively estimate the energy released during 
the giant burst to be: 

AXP IE 2259+586 E burst ~ 1.0 10 40 ergs . (4.30) 

On 1998 August 27, a giant burst from the soft gamma ray repeater SGR 1900+14 was 
recorded. The estimate energy released was: 

SGR 1900+14 E burst ~ 1.0 10 44 ergs . (4.31) 

As we have already discussed in Sect. the energy released during a burst in magnetars 
is given by the magnetic energy directly injected and dissipated into the magnetosphere, 
Eq. (ETra . We rewrite Eq. (ETTH1) as 

SE ext „ I o3 B 2 ^Bs_ _ 44 ( B S V 5Jh ( , 

oE B - s R B S Bs _ 2.6 10 ergs I 10 h Q auss J Bs ' { } 
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So that, combining Eqs. (|Q2l . (jOTj) . (PT3UJ) and P~2Tf) we get: 



iS00-/-^ ~ ^ ~ 0.70 10- 2 , 

IE 2259+586 — ^ ~ — - ~ 0.27 10~ 2 . 

Thus, according to Eq. (|4.20j) we may estimate the sudden variation of z>: 

^ 2 ^£ ~ 10- 2 , (4.34) 
v R c 

for both glitches. On the other hand we have: 



5I\ I5e c -e ( R c \ 5 5 R c 15 / R c \ 5 5 R c 



I J „„„„ 2 £ \ R j R r 2 \ R ) R r 

leading to: 



(4.35) 



SCR 1900+14 ( y ) ~ 3.31 10" 



i 



■ 



AXP IE 2159 i 586 f y ] ~ 1.31 10" 7 . 



(4.36) 



On the other hand, we expect that during the giant glitch {^surf ~ ^ 6 ' ^ s a con " 
sequence, for AXP IE 2259+586 the core contribution is negligible with respect to the 
surface contribution to 5L In other words, AXP IE 2259+586 displays a normal glitch 
with £ ~ 10~ 6 . On the contrary, Eq. (IQfij) indicates that SGR 1900+14 suffered a 
braking glitch with — ~ —3.34 10~ 4 giving: 



V 

1-3 



A P ~ 1.72 10~ d sec . (4.37) 

We would like to stress that our theory is in remarkable agreement with observations, for 
a glitch of size ^ = 4.24(11) 10~ 6 was observed in AXP IE 2259+586 which preceded 
the burst activity [HOI- Moreover, our theory predicts a sudden increase of the spin-down 
torque according to Eq. (|4.34|) . In Ref. [60] it is pointed out that it was not possible to give 
a reliable estimate of the variation of the frequency derivative since the pulse profile was 
undergoing large changes, thus compromising the phase alignment with the pulse profile 
template. Indeed, as discussed in Sect. 15.11 soon after the giant burst AXP IE 2259+586 
suffered an intense burst activity. Now, according to our theory, during the burst activ- 
ity there is both a continuous injection of magnetic energy into the magnetosphere and 
variation of the magnetic torque explaining the anomalous timing noise observed in IE 
2259+586. In addition, the authors of Ref. jHI] reported a gradual increase of the nomi- 
nal spin-down rate and a discontinuous spin down event associated with the 1998 August 
27 flare from SGR 1900+14- Extrapolating the long-term trends found before and after 
August 27, they found evidence of a braking glitch with AP ~ 0.57(2) 10~ 3 sec. In view 
of our theoretical uncertainties, the agreement with our Eq. (J4.37|) is rather good. 
We feel that it is worthwhile to point out that the standard magnetar theory is completely 
unable to predict the remarkable evidence of braking glitches. As a matter of fact, to our 
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knowledge, the only attempt to explain the braking glitch observed in SGR 1900+14 is 
done in Ref. jlSj where it is suggested that violent August 27 event involved a glitch. 
However, the magnitude of the glitch was estimated by scaling to the largest glitches in 
young, active pulsars with similar spin-down ages and internal temperature. In this way 
they deduced the estimate |^| ~ 1CT 5 to 10~ 4 . In our opinion, this can hardly be consid- 
ered a valid explanation for the braking glitch. First, the authors of Ref. overlooked 
the well known fact that radio pulsars display normal glitches and no braking glitches. 
Second, these authors cannot explain why AXP IE 2259+586 displayed a normal glitch 
instead of a braking glitch. 

Let us conclude this Section by briefly discussing the 2004 December 27 giant flare from 
SGR 1806-20. During this tremendous outburst SGR 1806-20 released a huge amount of 
energy, Eb urst ~ 10 46 ergs. Using the spin-down parameters reported in Ref. 

SGR 1806-20 P ~ 7.55 sec , P ~ 5.5 10~ 10 , B s ~ 2.0 10 15 Gauss , (4.38) 

we find: 

SGR 1806-20 ^ ~ ^ ~ 9.6 10^ 2 . (4.39) 

Thus, we predict that SGR 1806-20 should display a gigantic braking glitch with ^ ~ 
2.4 HT 2 , or : 

AP ~ 1.8 lO" 1 sec . (4.40) 



4.2 QUIESCENT LUMINOSITY 

The basic mechanism to explain the quiescent X-ray emission in our magnetars is the 
internal dissipation of magnetic energy. Our mechanism is basically the same as in the 
standard magnetar model based on neutron star 05] . Below we shall critically compare 
our proposal with the standard theory. 

In Section 0] we showed that during a glitch there is a huge amount of magnetic energy 
released into the magnetar: 

-6E%* ~ -R?B%(— V^£. (4.41) 
B 3 3 \Rj R c 1 ; 

As in previous Section, we use SGR 1900+14 and AXP IE 2259+586 as prototypes for 
soft gamma ray repeaters and anomalous X-ray pulsars, respectively. Using the results 
of Sect. 14.11 we find: 



SGR 1900+14 - 6 Eg* ~ 3.0 10 47 erg 
AXP IE 2259+586 - 5 E l £ l ~ 4.8 10 45 erg 



SRc 

Rc 
5R r 

R f 



(4.42) 



This release of magnetic energy is dissipated leading to observable surface luminosity. To 
see this, we need a thermal evolution model which calculates the interior temperature dis- 
tribution. In the case of neutron stars such a calculation has been performed in Ref. [HI] , 
where it is showed that the isothermal approximation is a rather good approximation in 
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the range of inner temperatures of interest. The equation which determines the thermal 
history of a p-star has been discussed in Ref. [H] in the isothermal approximation: 



Cy ^ = - (L v + L 7 ) , (4.43) 



where L v is the neutrino luminosity, L 7 is the photon luminosity and Cy is the spe- 
cific heat. Assuming blackbody photon emission from the surface at an effective surface 
temperature T$ we have: 

L 7 = 4nR 2 a SB T$, (AAA) 

where <jsb is the Stefan — Boltzmann constant. In Ref. |9 i we assumed that the surface 
and interior temperature were related by: 

^ = 10- 2 a . (4.45) 

Equation (|4.45j) is relevant for a p-star which is not bare, namely for p-stars which are 
endowed with a thin crust. It turns out that p-stars have a sharp edge of thickness 
of the order of about 1 fermi. On the other hand, electrons which are bound by the 
coulomb attraction, extend several hundred fermis beyond the edge. As a consequence, 
on the surface of the star there is a positively charged layer which is able to support a thin 
crust of ordinary matter. The vacuum gap between the core and the crust, which is of 
order of several hundred fermis, leads to a strong suppression of the surface temperature 
with respect to the core temperature. The precise relation between Ts and T could be 
obtained by a careful study of the crust and core thermal interaction. In any case, our 
phenomeno logical relation Eq. (|4.45|) allows a wide variation of Ts, which encompasses the 
neutron star relation (see, for instance, Ref. [03)- Moreover, our cooling curves display a 
rather weak dependence on the parameter a in Eq. (|4.45|) . Since we are interested in the 
quiescent luminosity, we do not need to known the precise value of this parameter. So 
that, in the following we shall assume a ~ 1. In other words, we assume: 

T s ~ 1(T 2 T . (4.46) 

Obviously, the parameter a is relevant to evaluate the surface temperature once the core 
temperature is given. Note that, in the relevant range of core temperature T ~ 10 8 °K, 
our Eq. (J4.46)) is practically identical to the parametrization adopted in Ref. jlH] within 
the standard magnetar model: 

T S * 1.3 10 6 °K (^) 9 • (4-47) 

The neutrino luminosity L u in Eq. ()4.43|) is given by the direct /5-decay quark reactions, 
the dominant cooling processes by neutrino emission. From the results of Ref. ||, we find: 

L v ~ 1.22 10 36 ^j- T 9 8 , (4.48) 

where T 9 is the temperature in units of 10 9 °K. Note that the neutrino luminosity L v 
has the same temperature dependence as the neutrino luminosity by modified URCA 
reactions in neutron stars (see, for instance Ref. [30 J, but it is more than two order of 
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c 
c 



magnitude smaller. From the cooling curves reported in Ref. [H] we infer that the surface 
and interior temperature are almost constant up to time r ~ 10 5 years. Observing that 
magnetars candidates are rather young pulsar with r age < 10 5 years, we may estimate the 
average surface luminosities as: 

L 7 ~ _^£§_ . (4.49) 

T age 

We assume r age w r c for SGR 1900+ 14- On the other hand, as discussed in Section [H we 
known that for AXP IE 2259+586 r age ~ 10 3 years < r c . We get: 

SGR 1900+14 L 7 ^ 1.3 1 37 ^ ^ 
.4XP IE 2259+586 L~ ~ 1.5 1 35 — — - . 

So that it is enough to assume that SGR 1900+14 suffered in the past a glitch with 
^ ~ 10~ 2 to sustain the observed luminosity L 7 ~ 10 35 ^j- (assuming a distance of 
about 10 kpc). In the case of AXP IE 2259+586, assuming a distance of about 3 kpc, 
the observed luminosity is L 7 ~ 10 34 — , so that we infer that this pulsar had suffered in 
the past a giant glitch with ^ ~ 10 _1 , quite similar to the recent SGR 1806-20 glitch. 
Let us discuss the range of validity of our approximation. Equation (J4.49|) is valid as 
long as L 7 dominates over L v , otherwise the star is efficiently cooled by neutrino emission 
and the surface luminosity saturates to L™ ax . We may quite easily evaluate this limiting 
luminosity from L™ ax ~ L u . Using Eq. ()4.46|) and R ~ 10 Km, we get: 

L raax _ 4 3 1Q 37 _ U ^ 

7 S 

Note that, since our neutrino luminosity is reduced by more than two order of magnitude 
with respect to neutron stars, L™ ax is about two order of magnitude greater than the 
maximum allowed surface luminosity in neutron stars [01] . Thus, while our theory allows 
to account for observed luminosities up to 10 36 — , the standard model based on neutron 
stars is in embarrassing contradiction with observations. 

Let us, finally, comment on the quiescent thermal spectrum in our theory. As already 
discussed, the origin of the quiescent emission is the huge release of magnetic energy 
in the interior of the magnetar. Our previous estimate of the quiescent luminosities 
assumed that the interior temperature distribution was uniform. However, due to the 
huge magnetic field, the thermal conductivity is enhanced along the magnetic field. This 
comes out to be the case for both electron and quarks, since we argued that the magnetic 
and chromomagnetic fields are aligned . As a consequence, we expect that the quiescent 
spectrum should be parameterized as two blackbodies with parameter R\ , T\ and R2 , T2, 
respectively. Since the blackbody luminosities L 7 i and L 7 2 are naturally of the same 
order, our previous estimates for the quiescent luminosities are unaffected. Moreover, 
since the thermal conductivity is enhanced along the magnetic field, the high temperature 
blackbody, with temperature T 2 , originates from the heated polar magnetic cups. Thus 
we have: 

Ri < R , R 2 < IKm , 

~ 4.52 
T 2 > 7i , R 1 > R 2 , 



20 



Note that there is a natural anticorrelation between blackbody radii and temperatures. 
Customary, the quiescent spectrum of anomalous X-ray pulsars and soft gamma ray 
repeaters is fitted in terms of blackbody plus power law. In particular, it is assumed 
that the power law component extends to energy greater than an arbitrary cutoff energy 
E C utoff — 2 KeV. It is worthwhile to stress that these parameterizations of the quiescent 
spectra are in essence phenomenological fits, for there are not sound physical motiva- 
tions. Indeed, within the standard magnetar model jlHj the power law should be related 
to hydromagnetic wind accelerated by Alfven waves. However, any physical justification 
for the arbitrary cutoff energy E cuto ff is lacking. Moreover, the luminosity of the wind 
emission should increase with magnetic field strength as L wind ~ L PL ~ Bg. On the other 
hand, the blackbody luminosities should scale as B S A [33]. So that the ratio L PL /L BB 
decreases with increasing magnetic field strengths, contrary to observations [66 . Finally, 
observations of a small energy dependence of pulsed fraction in some anomalous X-ray 
pulsars requires ad hoc tuning of the blackbody and power law components. Thus, we 
see that the standard magnetar model is in striking contradictions with observations. 
On the contrary, in our theory well defined physical arguments lead to the two black- 
body representation of the quiescent spectra, whose parameters are constrained by our 
Eq. ([4.52)1 . As a matter of fact, we have checked in literature that the quiescent spectra of 
both anomalous X-ray pulsars and soft gamma ray repeaters could be accounted for by 
two balckbodies. For instance, in Ref. jHZI the quiescent spectrum of AXP IE 1841-045 
is well fitted with the standard power law plus blackbody (reduced x 2 /dof ~ 1-11), nev- 
ertheless the two blackbody model gives also a rather good fit (reduced y 2 /dof ~ 1.12). 
Interestingly enough the blackbody parameters: 

R 1 = 5.7 ±S:f Km , 7\ = 0.47 ± 0.02 KeV , 

R 2 = 0.3618$ Km > T 2 = 1.5 la? KeV , ( " >> 

are in agreement with Eq. ([4.52)1 . Moreover, assuming that the power law component in 
the standard parametrization of quiescent spectra account for the hot blackbody com- 
ponent in our parametrization, we find that the suggestion L 7l ~ L 72 in Eq. ([4.52)1 is 
in agreement with observations [SH]- It should be stressed, however, that the two black- 
bodies are not the whole story. In Appendix we show that thermal photons originating 
from the hot polar cups are reprocessed by electrons trapped above the polar cups. These 
electrons, which are responsible for the faint low energy spectrum, could result in broad 
spectral features in the quiescent spectrum. These spectral features, in turn, could result 
in observable deviations from the two blackbody fit. 

Another interesting consequence of the anisotropic distribution of the surface tempera- 
ture due to strong magnetic fields is that the thermal surface blackbody radiation will 
be modulated by the stellar rotation. As a matter of fact, in Ref. (HE! it is argued that 
the observed properties of anomalous X-ray pulsars can be accounted for by magnetars 
with a single hot region. It is remarkable that our interpretation explains naturally the 
observed change in pulse profile of SGR 1900+14 following the 1998 August 27 giant 
flare. In addition, the thermal radiation reprocessed by electrons near the polar cups 
could result in an effective description with two hot spots. It seem that our picture is in 
fair qualitative agreement with several observations. However, any further discussion of 
this matter goes beyond the aim of the present paper. 
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4.3 BURSTS 



In the present Section we discuss how glitches in our magnetars give rise to the burst 
activity from soft gamma-ray repeaters and anomalous X-ray pulsars. We said in Sect. H] 
that the energy released during a burst in a magnetar is given by the magnetic energy 
directly injected into the magnetosphere, Eq. (j4.19J) . Before addressing the problem of 
the dissipation of this magnetic energy in the magnetosphere, let us discuss what are the 
observational signatures at the onset of the burst. Observations indicate that at the onset 
of giant bursts the flux displays a spike with a very short rise time ti followed by a rapid 
but more gradual decay time t 2 . According to our previous discussion, the onset of bursts 
is due to the positive variation of the surface magnetic field 5B$, which in turn implies 
an sudden increase of the magnetic energy stored in the magnetosphere. Now, according 
to Eq. ()4.18j) we see that almost all the magnetic energy is stored in the region: 

R < r < 10 R . (4.54) 

So that the rise time is essentially the time needed to propagate in the magnetosphere 
the information that the surface magnetic field is varied. So that we are lead to: 

ti ~ 9 R ~ 3 1(T 4 sec , (4.55) 

which indeed is in agreement with observations. On the other hand, in our proposal the 
decay time t 2 depends on the physical properties of the magnetosphere. Indeed, it is 
natural to identify t 2 with the time needed to the system to react to the sudden variation 
of the magnetic field. In other words, we may consider the magnetosphere as a huge 
electric circuit which is subject to a sudden increase of power from some external device. 
The electric circuit reacts to the external injection of energy within a transient time. So 
that, in our case the time t 2 is a function of the geometry and the conducting properties 
of the magnetosphere. In general, it is natural to expect that t\ <C t 2 so that the time 
extension of the initial spike is: 

S t S pi ke ~ t 2 - h ~ t 2 ■ (4.56) 

Remarkably, observations shows that the observed giant bursts are characterized by almost 
the same 5 t spi ke- 

$ tspike ^ h - 0.1 sec , (4.57) 

signalling that the structure of the magnetosphere of soft gamma-ray repeaters and 
anomalous X-ray pulsars are very similar. Since the magnetic field is varied by SBs 
in a time St sp ik e , then from Maxwell equations it follows that it must be an induced 
electric field. To see this, let us consider the dipolar magnetic field in polar coordinate: 

2 B? R 3 cos 9 



B r 



B s R 3 sinO (4.58) 
Bq = 



B in = 



Thus, observing that is the time derivative of the magnetic field it is easy to find 



the induced azimuthal electric field: 



5B S R 3 sm9 „ , A . 

E v = +—± — , r > R . (4.59) 

spike ' 
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To discuss the physical effects of the induced azimuthal electric field Eq. ()4.59|h it is 
convenient to work in the co-rotating frame of the star. We assume that the magneto- 
sphere contains a neutral plasma. Thus, we see that charges are suddenly accelerated by 
the huge induced azimuthal electric field and thereby acquire an azimuthal velocity 
v v ~ 1 which is directed along the electric field for positive charges and in the opposite 
direction for negative charges. Now, it is well known that relativistic charged particles 
moving in the magnetic field B (r ) , Eq. (j4.58)) , will emit synchrotron radiation (HHI • As we 
discuss below, these processes are able to completely dissipate the whole magnetic energy 
injected into the magnetosphere following a glitch. However, before discuss this last point 
in details, we would like to point out some general consequences which lead to well defined 
observational features. As we said before, charges are accelerated by the electric field E v 
thereby acquiring a relativistic azimuthal velocity. As a consequence, they are subject to 
the drift Lorentz force F = qv^ x B(r), whose radial component is: 

3 



F r = -qv^Be ~ +qv v B s sin 6 ( - ) , (4.60) 

while the 9 component is: 

F e = +qv v B r ~ -2qv v B s cose(j \ . (4.61) 

The radial component F r pushes both positive and negative charges radially outward. 
Then, we see that the plasma in the outermost part of the magnetosphere is subject to 
a intense radial centrifugal force, so that the plasma must flow radially outward giving 
rise to a blast wave. On the other hand, Fg is centripetal in the upper hemisphere and 
centrifugal in the lower hemisphere. As a consequence, in the lower hemisphere charges 
are pushed towards the magnetic equatorial plane cos# = 0, while in the upper hemi- 
sphere (the north magnetic pole) the centripetal force gives rise to a rather narrow jet 
along the magnetic axis. As a consequence, at the onset of the giant burst there is an 
almost spherically symmetric outflow from the pulsar together with a collimated jet from 
the north magnetic pole. Interestingly enough, a fading radio source has been seen from 
SGR 1900+14 following the August 27 1998 giant flare jTUj. Indeed, the radio afterglow 
is consistent with an outflow expanding subrelativistically into the surrounding medium. 
This is in agreement with our model once one takes into account that the azimuthal elec- 
tric field is rapidly decreasing with the distance from the star, so that v v < 1 for the 
plasma in the outer region of the magnetosphere. However, we believe that the most 
compelling evidence in favour of our proposal comes from the detected radio afterglow 
following the 27 December 2004 gigantic flare from SGR 1806-20 HQ E21 E21 El E5]. 
Indeed, the fading radio source from SGR 1806-20 has similar properties as that observed 
from SGR 1900+14, but much higher energy. The interesting aspect is that in this case 
the spectra of the radio afterglow showed clearly the presence of the expected spherical 
non relativistic expansion together with a sideways expansion of a jetted explosion (see 
Fig. 1 of Ref. jTI] and Fig. 1 of Ref. [12]), i n striking agreement with our theory. Note 
that the standard magnetar model is unable to account for these observed features of the 
radio afterglow. Finally, the lower limit of the outflow E > 10 44 ' 5 ergs {75.. implies that 
the blast wave and the jet dissipate only a small fraction of the burst energy which is 
about 10 46 ergs (see Section f4.1|) . Thus, we infer that almost all the burst energy must 
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be dissipated into the magnetosphere. In the co-rotating frame of the star the plasma, 
at rest before the onset of the burst, is suddenly accelerated by the induced electric field 
thereby acquiring an azimuthal velocity v v ~ 1. Now, relativistic charges are moving in 
the dipolar magnetic field of the pulsar. So that, they will lose energy by emitting syn- 
chrotron radiation until they come at rest. Of course, this process, which involves charges 
that are distributed in the whole magnetosphere, will last for a time tdi S much longer 
that St S pike- Actually, will depend on the injected energy, the plasma distribution and 
the magnetic field strength. Moreover, one should also take care of repeated charge and 
photon scattering. So that it is not easy to estimate t dis without a precise knowledge of 
the pulsar magnetosphere. At the same time, the fading of the luminosity with time, the 
so-called light curve L(t), cannot be determined without a precise knowledge of the mi- 
croscopic dissipation mechanisms. However, since the dissipation of the magnetic energy 
involves the whole magnetosphere, it turns out that we may accurately reproduce the 
time variation of L(t) without a precise knowledge of the microscopic dissipative mecha- 
nisms. Indeed, in Sect. El we develop an effective description where our ignorance on the 
microscopic dissipative processes is encoded in few macroscopic parameters, which allows 
us to determine the light curves. In the remaining of the present Section we investigate 
the spectral properties of the luminosity. To this end, we need to consider the synchrotron 
radiation spectral distribution. Since radiation from electrons is far more important than 
from protons, in the following we shall focus on electrons. It is well known that the 
synchrotron radiation will be mainly at the frequency [TE] (see also Ref. 



2 — , 4.62 

rrif, 



where 7 is the electron Lorentz factor. Using Eq. (|4.58|) we get: 

uj(r) ~ 7 2 — (-) , R < r < 10 R . (4.63) 
m e \rj 

It is useful to numerically estimate the involved frequencies. To this end, we consider the 
giant flare of 1998, August 27 from SGR 1900+14: 

B s ~ 7 .4 10 14 Gauss , ~ 1(T 2 . (4.64) 

Bs 

So that, from Eq. (|4.63|) it follows: 

u{r) ~ 7 2 8.67 MeV ( , R < r < 10 R , (4.65) 



or 

wi ~ 7 2 8.67 if el/ < u < lu 2 ~ 7 2 8.67 MeV . (4.66) 

The power injected into the magnetosphere is supplied by the azimuthal electric field 
during the initial hard spike. So that to estimate the total power we need to evaluate the 
power supplied by the azimuthal electric field. Let us consider the infinitesimal volume 
dV = r 2 sin OdrdOdip; the power supplied by the induced electric field E v in dV: 

dW Ev ~ n e e v v R 3 sin 2 6drd6d<p , (4.67) 

Tjspike 
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where n e is the electron number density. Since the magnetosphere is axially symmetric 
it follows that n e cannot depend on ip. Moreover, within our theoretical uncertainties we 
may neglect the dependence on 9. So that, integrating over 9 and <p we get: 

dW Eip ~ 2 it 2 n e e v v R 3 dr . (4.68) 



spike 

In order to determine the spectral distribution of the supplied power, we note that to a 
good approximation all the synchrotron radiation is emitted at u m , Eq. ()4.62|) . So that, 
we may use Eq. (I4.65j) to get: 

dr ~ tL 7 § ( £^ ) _L du . (4.69) 

UJ 3 



3 \ m 



Inserting Eq. (|4.69|) into Eq. (|4.68|) we obtain the spectral power: 

/ x , 2tt 2 5£ s 4 2 feB s \ 1/3 1 , . , 

^ ~ — n e e — — v v R* 73 — i — rfcu , (4.70) 

3 Olgpike V ?7le / o;3 

while the total luminosity is given by: 

/•U>2 

L — F(u) du . (4.71) 



Note that L is the total luminosity injected into the magnetosphere during the initial 
hard spike. So that, since the spike lasts 5 t sp ik e , we have: 



J burst 



8 t spike L , (4.72) 



where Eb urst is the total burst energy. In the case of the 1998 August 27 giant burst from 
SGR 1900+14 the burst energy is given by Eq. (JOTJ). Thus, using Eqs. (ET721 and (ED77f) 

we have: 

L ~ 10 45 ^ , (4.73) 

sec 

which, indeed, is in agreement with observations. It is worthwhile to estimate the electron 
number density needed to dissipate the magnetic energy injected in the magnetosphere. 
To this end, we assume an uniform number density. Thus, using Eqs. (14.71)1 . (|4.70J) and 
(JHIEJ) we get: 

L ~ 18 vr 2 n e e R 4 , (4.74) 

spike 

where we used v v ~ 1. Specializing to the August 27 giant burst we find: 

n e ~ 2.0 10 14 cm" 3 , (4.75) 

indeed quite a reasonable value. Soon after the initial spike, the induced azimuthal 
electric field vanishes and the luminosity decreases due to dissipative processes in the 
magnetosphere. As thoroughly discussed in Sect. El it is remarkable that the fading of the 
luminosity can be accurately reproduced without a precise knowledge of the microscopic 
dissipative mechanisms. So that combining the time evolution of the luminosity L(t), 
discussed in Sect. El with the spectral decomposition we may obtain the time evolution 
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of the spectral components. In particular, firstly we show that starting from Eq. ()4.70j) 
the spectral luminosities can be accounted for by two blackbodies and a power law. After 
that, we discuss the time evolution of the three different spectral components. 
The spectral decomposition Eq. ()4.70j) seems to indicate that the synchrotron radiation 
follows a power law distribution. However, one should take care of reprocessing effects 
which redistribute the spectral distribution. To see this, we note that photons with energy 
uj > 2m e quickly will produce copiously almost relativistic e pairs. Now, following 
Ref. JUj, even if the particles are injected steadily in a time St sp ik e , it is easy to argue that 
the energy of relativistic particles is rapidly converted due to comptonization to thermal 
photon-pair plasma. Since the pair production is quite close to the stellar surface, we may 
adopt the rather crude approximation of an uniform magnetic filed B ~ B$ throughout 
the volume V v i asraa ~ 12 7ri? 3 . Since typical magnetic fields in magnetars are well above 
Bqed-i electrons and positrons sit in the lowest Landau levels. In this approximation we 
deal with an almost one dimensional pair plasma whose energy density is 
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— — if Tpi asma \ 2 TTlf 



u e ~ m e (n e + +n e -) ~ - — -j eB s m e exp(-- ) , (4.76) 

(27t) 2 \ m e J l p i 

OS OKI 

for T p i asma < m e , Tpi asma being the plasma temperature. Thus, the total energy density 
of the thermal photon-pair plasma is: 

7T 2 

u = u e + m 7 ~ u e + — T 4 lasma . (4.77) 

The plasma temperature is determined by equaling the thermal energy Eq. (|4.77J1 with 
the fraction of burst energy released in the spectral region to > 2 m e . It is easy to find: 

Epa irs ~ 0.147 E burst , (4.78) 

where for the numerical estimate we approximated uii ~ lOKeV and 002 — lOMeV, 
corresponding to mildly relativistic electrons in the magnetosphere. So that we have: 

r> 2 fTplasma\^ , m e \ 4 Ep a i rs _„n 

™ e m eX P(~^ ) + Tr T P lasma ~ 77 ■ ( 4 -79) 

\ ifl/e J -L plasma * plasma 



(2vr)i 

In the case of August 27 giant burst from SGR 1900+14 we find: 

exp(-i) + 0.311 x 4 ~ 1.32 10" 2 , x = , (4.80) 

x m e 

whose solution gives Tpi aS ma — 135 KeV. However, this is not the end of the whole 
story. Indeed, our thermal photon-pair plasma at temperature T p i a sma will be reprocessed 
by thermal electrons on the surface which are at temperature of the thermal quiescent 
emission Tq < 1 KeV. So that, photons at temperature T pZasma ^> Tq are rapidly cooled 
by Thompson scattering off electrons in the stellar atmosphere, which extends over several 
hundreds fermis beyond the edge of the star. The rate of change of the radiation energy 
density is given by [77] : 

1 du~ 4 Ut no /m m , . , „ 

~TT~ — (T Q - Tp lasma ) , 4.81 

w-v at m P 
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where tiq is the number density of electrons in the stellar atmosphere. The electron 
number density in the atmosphere of a p-star is of the same order as in strange stars, where 
uq ~ 10 33 cm~ 3 (see for instance Ref [7Hj). So that, due to the very high electron density 
of electrons near the surface of the star, the thermal photon-pair plasma is efficiently 
cooled to a final temperature much smaller than T p i asma . At the same time, the energy 
transferred to the stellar surface leads to an increase of the effective quiescent temperature. 
Therefore we are lead to conclude that during the burst activity the quiescent luminosity 
must increase. Let T\ be the final plasma temperature, then we see that the thermal 
photon-pair plasma contribution to the luminosity can be accounted for with an effective 
blackbody with temperature T\ and radius R\ of the order of the stellar radius. As a 
consequence the resulting blackbody luminosity is: 

L x = ATtRla SB Tt , Ri < R . (4.82) 

In general, the estimate of the effective blackbody temperature T\ is quite difficult. How- 
ever, according to Eq. (|4.78|) we known that L\ must account for about 0.147 of the total 
luminosity. So that we have: 

L x {t) ~ 0.147 L(t) . (4.83) 

This last equation allow us to determine the blackbody temperature. For instance, soon 
after the hard spike we have L(0) ~ f. burst ~ 10 45 for the giant burst from SGR 

y ' utspike Sec ° 

1900+14. Thus, using R x ~ R, from Eq. we get: 

71(0) ~ 61 KeV , (4.84) 
with surface luminosity Li(0) ~ 10 44 ^Sl, 

Let us consider the remaining spectral power with uj < 2m e . We recall that the spectral 
power Eq. ()4.70|) originates from the power supplied by the induced electric field Eq. ()4.68|) . 
It is evident from Eq. (|4.68|) that, as long as v v ~ 1, the power supplied by the electric 
field Ep does not depend on the mass of accelerated charges. Since the plasma in the 
magnetosphere is neutral, it follows that protons acquire the same energy as electrons. 
On the other hand, since the protons synchrotron frequencies are reduced by a factor 
— , the protons will emit synchrotron radiation near u\. As a consequence, photons with 
frequencies near u\ suffer resonant synchrotron scattering, which considerably redistribute 
the available energy over active modes. On the other hand, for u ^> LO\ the spectral power 
will follows the power law Eq. (J4.70|) . Thus, we may write: 

F(u) ~ -4 , 5wi < w < 2m e , (4.85) 

where we have somewhat arbitrarily assumed the low energy cutoff ~ 5 u)\ . On the other 
hand, for uj < 5 U\ the redistribution of the energy by resonant synchrotron scattering 
over electron and proton modes lead to an effective description of the relevant luminosity 
as thermal blackbody with effective temperature and radius T 2 and R2, respectively. 
Obviously, the blackbody radius R2 is fixed by the geometrical constrain that the radiation 
is emitted in the magnetosphere at distances r < 10 R. So that we have: 

R 2 < 10 R . (4.86) 
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The effective blackbody temperature T 2 can be estimate by observing that the integral of 
the spectral power up to 5 uj\ account for about the 60 % of the total luminosity. Thus, 
we have: 

L 2 (t) ~ 0.60 L(t) , (4.87) 

where 

L 2 = 4nB*(T SB T* , R 2 < 10R . (4.88) 

Equations ()4.87j) and (|4.88j) can be used to to determine the effective blackbody temper- 
ature. If we consider again the giant burst from SGR 1900+14, soon after the hard spike, 
assuming R 2 — 10 R, we readily obtain: 

T 2 (0) ~ 27 KeV . (4.89) 

To summarize, we have found that the spectral luminosities can be accounted for by two 
blackbodies and a power law. In particular for the blackbody components we have: 

R x < R , R 2 < 10 R ; T 2 < T x , R x < R 2 

L x ~ 0.15 L , L 2 ~ 0.60 L . [ ' ' 

Interestingly enough, Eq. ()4.90|) displays an anticorrelation between blackbody radii and 
temperatures, in fair agreement with observations. Moreover, the remaining 25% of the 
total luminosity is accounted for by a power law leading to the high energy tail of the 
spectral flux: 

~j~ ~ E~ a , a ~ 2.33 , (4.91) 
dE 

extending up to E ~ 2m e ~ 1 MeV. Indeed, the high energy power law tail is clearly 
displayed in the giant flare from SGR 1900+14 (see Fig. 3 in Ref. [10]), and in the recent 
gigantic flare from SGR 1806-20 (see Fig. 4 in Ref. (HO! ) - 

It is customary to fit the spectra with the sum of a power law and an optically thin thermal 
bremsstrahlung. It should be stressed that the optically thin thermal bremsstrahlung 
model is purely phenomenological and without a physical basis. In view of this, a direct 
comparison of our proposal with data is problematic. Fortunately, the authors of Ref. [HI] 
tested several spectral functions to the observed spectrum in the afterglow of the giant 
outburst from SGR 1900+14- In particular they found that, in the time interval 68 sec < 
t < 195 sec, the minimum \ 2 spectral model were composed by two blackbody laws plus 
a power law. By fitting the time averaged spectra they reported 



T 2 ~ 9.3 KeV , T x ~ 20.3 KeV , a ~ 2.8 . (4.92) 

Moreover, it turns out that the power law accounts for approximately 10% of the total en- 
ergy above 25 KeV , while the low temperature blackbody component accounts for about 
85% of the total energy above 25 KeV . In view of our neglecting the contribution to 
energy from protons, we see that our proposal is in accordance with the observed energy 
balance. Unfortunately, in Ref. [HI] the blackbody radii are not reported. To compare our 
estimate of the blackbody temperatures with the fitted values in Eq. (14.92}) . we note that 
our values reported in Eqs. ()4.84j) and ()4.89j) correspond to the blackbody temperatures 
soon after the initial hard spike. Thus, we need to determine how the blackbody temper- 
atures evolve with time. To this end, we already argued that soon after the initial spike 
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the luminosity decreases due to dissipative processes in the magnetosphere. In Sect. El we 
show that the fading of the luminosity can be accurately reproduced without a precise 
knowledge of the microscopic dissipative mechanisms. In particular, the relevant light 
curve is given by Eqs. ()5.6|) and (j5.10|) . At t = we have seen that the total luminosity is 
well described by three different spectral components. During the fading of the luminos- 
ity, it could happens that microscopic dissipative processes modify the different spectral 
components. However, it is easy to argue that this does not happens. The crucial point is 
that the three spectral components originate from emission by a macroscopic part of the 
magnetosphere; moreover the time needed to modify a large volume of magnetosphere by 
microscopic processes is much larger than the dissipation time t<n s ~ 10 2 sec. Then we 
conclude that, even during the fading of the luminosity, the decomposition of the lumi- 
nosity into three different spectral components retain its validity. Now, using the results 
in Sect. 03 we find: 

L(t " 6 f g6c) * 3.67 10- , L(t " 1%gec) ^ 1.67 10- . (4.93) 
L(0) ' L(0) v ; 

Combining Eq. flOS]) with Eqs. (jPSjl . lQ7l) and KM we obtain: 



T 2 (t ~ 68sec) ~ 11.8 KeV , T 2 (t ~ 195 sec) ~ 9.7 KeV 
Tx(t ~ 68sec) ~ 26.7 KeV , T^t ~ 195sec) ~ 21.9 KeV 



(4.94) 

in reasonable agreement with Eq. ()4.92j) . Finally, let us comment on the time evolution of 
the spectral exponent a in the power law Eq. (|4.91j) . From Eq. 1)4.70)1 it follows that high 
energy modes have less energy to dissipate. Accordingly, once a finite amount of energy 
is stored into the magnetosphere, the modes with higher energy become inactive before 
the lower energy modes. As a consequence, the effective spectral exponent will increases 
with time and the high energy tail of the emission spectrum becomes softer, in perfect 
agreement with observations. This explains also why the fitted spectral exponent a in 
Eq. (j4.92j) is slightly higher than our estimate in Eq. (|4.91j) . 

4.4 HARDNESS RATIO 

Recently, it has been reported evidence for a hardness-intensity anti correlation within 
bursts from SGR 1806-20 [H2|- Indeed, the authors of Ref. jH2] reported observations of 
the soft gamma ray repeaters SGR 1806-20 obtained in October 2003, during a period of 
bursting activity. They found that some bursts showed a significant spectral evolution. 
However, in the present Section we focus on the remarkable correlation between hardness 
ratio and count rate. Following Ref. jH2j we define the hardness ratio as: 

M = ft! ■ (495) 

where H and S are the background subtracted counts in the ranges 40 — 100 KeV and 
15 — 40 KeV respectively. In Figure 0] we report the hardness ratio data extracted from 
Fig. 3 of Ref. [B2]- A few comments are in order. First, the hardness ratio becomes 
negative for large enough burst intensities. Moreover, there is a clear decrease of the 
hardness ratio with increasing burst intensities. Note that, no detailed predictions are 
available within the standard magnetar model. On the other hand, within our approach 
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Figure 4: Hardness- intensity plot of the time resolved hardness ratio, Eq. (|4.95jl . Data 
have been extract from Fig. 3 of Ref. |H21- Blue lines are our Eqs. (|4.1UU|) and (j4.1(J2j) . 



we are able to explain why the hardness ratio is negative and decreases with increasing 
burst intensities. To see this, we note that the hardness ratio Eq. ()4.95|) is defined in 
terms of total luminosities in the relevant spectral intervals. Thus, to determine the total 
luminosity in the spectral interval oj\ — uj 2 we may use: 

L{u x -u 2 ) = / F(u) duo , (4.96) 

where F(u) is given by Eq. (J4.7Uj) . A straightforward integration gives: 

2 5Bs d4 / r i r 2 



L{lo\ — lu 2 ) ~ 2it I n e e— — [- - -) , (4.97) 

Ot spike V -K ti' 



where r\ and r 2 are given by: 



2 eB s ( R 



wu ~ 7 — • (4-98) 

m e V r i,2, 



Assuming 7 ~ 1, we may rewrite Eq. 1)4.98)1 as: 

u x ,i ^ 10 MeV ■ ( 4 - 99 ) 

Using Eqs. ()4.97j) and ()4.99j) it is easy to determine the hardness ratio: 

m-lOOKeV) - L(15-40KeV) 1.66 - 2.44 _ 

L(40 - 100 KeV) + L(15 - 40 KeV) 1.66 + 2.44 ' v ' ; 
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In Figure E| we display our estimate of the hardness ratio Eq. (j4.1UU|) . We see that data 
are in quite good agreement with Eq. (|4.1UU|) at least up to count rate ~ 5 10 3 counts . For 
larger count rates data seem to lie below our value. We believe that, within our approach, 
there is a natural explanation for this effect. Indeed, for increasing count rates we expect 
that the hard tail uj > 2 m e ~ 1 MeV of the spectrum will begin to contribute to the 
luminosity. According to the discussion in Sect. 14.31 these hard photons are reprocessed 
leading to an effective blackbody with temperature T%. Now, for small and intermediate 
bursts the blackbody temperature T\ is considerably smaller than Eq. ()4.84j) . so that the 
effective blackbody contributes mainly to the soft tail of the spectrum. Obviously, the 
total luminosity of the effective blackbody is: 

5B 

L(l-10MeV) ~ 2n 2 n e e- — — i? 4 (2.15 - 1) . (4.101) 

T spike 

Since this luminosity contributes to the soft part of the emission spectrum, Eq. ()4.1U0j) 
gets modified as: 

HR ~ L66 ~ 3,59 ~ - 0.37 . (4.102) 
1.66 + 3.59 v ' 

Equation (14.10211 is displayed in Fig. H for rates > 5 10 3 Note that we did not 

take into account the proton contribution to the luminosity. Observing that protons 
contribute mainly to luminosities at low energy uj < 10 KeV, we see that adding the 
proton contributions leads to smaller hardness ratios bringing our estimates to a better 
agreement with data. In any case, we see that our theory allows to explain in a natural 
way the puzzling anti correlation between hardness ratio and intensity. 

5 LIGHT CURVES 

In our magnetar theory the observed burst activities are triggered by glitches which inject 
magnetic energy into the magnetosphere where, as discussed in previous Sections, it is 
dissipated. As a consequence the observed luminosity is time depended. In this Section 
we set up an effective description which allows us to determine the light curves, i.e. the 
time dependence of the luminosity. In general, the energy injected into the magnetosphere 
after the glitch decreases due to dissipative effects described in Sect. 14 .3| leading to the 
luminosity L{t) = — dE ^ . Actually, the precise behavior of L{t) is determined once the 
dissipation mechanisms are known. However, since the dissipation of the magnetic energy 
involves the whole magnetosphere, we may accurately reproduce the time variation of 
L(t) without a precise knowledge of the microscopic dissipative mechanisms. Indeed, on 
general grounds we have that the dissipated energy is given by: 

L(t) = - = «(*) & , n < i , (5.1) 

where i] is the efficiency coefficient. Obviously the parameter n(t) does depend on the 
physical parameters of the magnetosphere. For an ideal system, where the initial injected 
energy is huge, the linear regime, where rj = 1, is appropiate. Moreover, we expect that the 
dissipation time ~ - is much smaller than the characteristic time needed to macroscopic 
modifications of the magnetosphere. Thus, we may safely assume n(t) ~ n constant. So 
that we get: 

L(t) = - ^ ~ n E . (5.2) 
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It is straightforward to solve Eq. (|5.2|) : 

E{t) = E exp(--) , L(t) = L exp(--) , L = — , r = — . (5.3) 

TO T T K 

Note that the dissipation time r = — encodes all the physical information on the mi- 
croscopic dissipative phenomena. Since the injected energy is finite, the dissipation of 
energy degrades with the decreasing of the available energy. Thus, the relevant equation 
is Eq. (|5.1|) with 77 < 1. In this case, solving Eq. (J5.1)) we find: 

V 

t \ 



L(t) = L 1 - — , (5.4) 

V t dis J 

where we have introduced the dissipation time: 

1 E^ T ' 

tdis = — • (5.5) 

Ko 1-1] 

Then, we see that the time evolution of the luminosity is linear up to some time tb rea k, 
after that we have a break from the linear regime 77 = 1 to a non linear regime with 77 < 1 . 
If we indicate with the total dissipation time, we get: 



L(t) = L exp( -) , < t < t break , 

TO 

L(t) = L{t U reak) I 1 — - - I , tbreak <■ t < tdis 

\ ^dis tbreak J 



(5.6) 



Equation (|5.6jl is relevant to describe the light curve after a giant burst, where there is a 
huge amount of magnetic energy dissipated into the magnetosphere. It is interesting to 
compare our light curves, Eq. (|5.6j) . with the standard magnetar model. The decay of the 
luminosity in the standard magnetar model is due to the evaporation by a fireball formed 
after a giant burst and trapped onto the stellar surface jlHl EH] • Indeed, the authors of 
Ref. [SI] considered the light curves after the giant flare of 1998, August 27 from SGR 
1900+14, and the giant flare of 1979 March 5 from SGR 0526-66. Assuming that the 
luminosity varies as a power of the remaining fireball energy L ~ E a , they found: 



l-a 



L(t) = L 1 - — , (5.7) 

\ ^evap J 

where t evap is the time at which the fireball evaporates, and the index a accounts for the 
geometry and the temperature distribution of the trapped fireball. For a spherical fireball 
of uniform temperature a = |, so that the index a must satisfies the constrain: 

2 , s 
a ^ - . 5.8 

3 K J 

Note that our Eq. ()5.6|) reduces to Eq. ()5.7|) if t orea k = and r\ = a. However, we stress 
that our efficiency exponent must satisfy the milder constraint 77 ^ 1. 

The authors of Ref. [HI] performed a best fit of the light curve of the August 27 flare, 
background subtracted and binned to 5 sec, to Eq. (J5.7|) and found: 

a = 0.756 ± 0.003 , t evap = 414 sec . (5.9) 
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Figure 5: Light curve after the giant flare of 1998, August 27 from SGR 1900+14- 
Red continuous line is the light curve in the standard magnetar model, Eq. ()5.7|) with 
parameters given in Eq. (|5.9|) . Blue line is our light curve Eq. ()5.6|) with parameters in 

Eq. flszrap 

Indeed, from Fig. 2 in Ref. [SJ one sees that the trapped fireball light curve account for 
the decay trend of the experimental light curve and matches the sudden final drop of the 
flux. However, it should stressed that the fitting parameter a in Eq. ()5.9|) does not satisfy 
the physical constraint Eq. (|5.8|) . Even more, any deviations from spherical geometry or 
uniform temperature distribution lead to parameters a smaller than the upper bound |. 
Moreover, the trapped fireball light curve underestimate by about an order of magnitude 
the measured flux during the first stage of the outburst. We interpreted the different 
behavior of the flux during the initial phase of the outburst as a clear indication of the 
linear regime described by our Eq. (|5.3|) . As a matter of fact, we find that the measured 
light curve could be better described by Eq. (|5.6|) with parameters (see Fig. |5J): 

t = 8.80 sec , U rea k = 20 sec , r] = 0.756 , tdis = 414 sec . (5.10) 

The same criticisms apply to the fit within the standard magnetar model of the light 
curve after the giant flare of 1979 March 5 from SGR 0526-66. The trapped fireball light 
curve fit in Ref. jHI] gives: 

a = 0.71 ± 0.01 , t evap = 163 ± 5 sec . (5.11) 

Again the parameter a exceeds the bound Eq. (j5.8J) . and the fit underestimates the flux 
during the first stage of the outburst (see Fig. 14, Ref. [EI])- Fitting our Eq. ()5.6|) to the 
measured flux reported in Fig. 14 of Ref. [HI], we estimate: 

r = 15 sec , tb rea k = 20 sec , t] = 0.71 , tdi s = 163 sec . (5-12) 
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Figure 6: Light curve after the giant flare of 1979 March 5 from SGR 0526-66. Red con- 
tinuous line is the light curve in the standard magnetar model, Eq. ()5.7|) with parameters 
given in Eq. ([5.110 . Blue line is our light curve Eq. (|5.6[) with parameters in Eq. (J5.12[) 



In Figs. and El we compare our light curves Eqs. ([5.60 . (J5.1U|) and ([5.120 with the best 
fits performed in Ref. [HI]. Obviously, both light curves agree for t > U rea k, while in the 
linear regime t < Ureak, where the trapped fireball light curves underestimate the flux, 
our light curve follow the exponential decay and seem to be in closer agreement with 
observational data. 

Several observations indicate that after a giant burst there are smaller and recurrent 
bursts. According to our theory these small and recurrent bursts are the effect of several 
small glitches following the giant glitch. We may think about these small bursts like the 
seismic activity following a giant earthquake. These seismic glitches are characterized by 
light curves very different from the giant burst light curves. In the standard magnetar 
model these light curves are accounted for with an approximate t _ar decay [H3j. Below we 
compare this prediction with observations and argue that the standard theory is unable 
to adequately describe observational data. On the other hand, within our theory there is 
a natural way to describe the seismic burst activity. Indeed, during these seismic bursts, 
that we shall call settling bursts, there is an almost continuous injection of energy into the 
magnetosphere which tends to sustain an almost constant luminosity. This corresponds 
to an effective k in Eq. (|5.1[) which decreases smoothly with time. The simplest choice is: 

«(«) = ^_ . ( 5 . 13 ) 



34 



Inserting into Eq. (j5.1J) and integrating, we get: 



E(t) 



So that the luminosity is: 



-(1-V)— ln(l + M) 

K X 



1-77 



L(t) 



U 



(l + «it)" 
After defining the dissipation time: 



1 - (1 - v) — ln(l + Kit) 



i-n 



ln(l + KiUi 



K\ E 



i-v 



we rewrite Eq. (|5.15jl as 

L(t) 



[1 + Kit)' 



1 - 



Kq 1 - 7] 



ln(l + Kit) 
ln(l + Kit dis 



v 
1-1 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



Note that the light curve Eq. f)5 . 1 7j) depends on two characteristic time constants ^- and 
tdis- We see that Kitdis, which is roughly the number of small bursts occurred in the 
given event, gives an estimation of the seismic burst intensity. Moreover, since during the 
seismic bursts the injected energy is much smaller than in giant bursts, we expect that 
fitting Eq. (J5.17j) to the observed light curves will result in values of rj smaller than the 
typical values in giant bursts. In the following Sections we show that, indeed, our light 
curves Eq. ()5.17|) are in good agreement with several observations. 



5.1 AXP IE 2259+586 

On 2002, June 18 SGR-like bursts was recorded from AXP IE 2259+586. Coincident with 
the burst activity were gross changes in the pulsed flux, persistent flux, energy spectrum, 
pulse profile and spin down of the source |S0j . As discussed in previous Sections, these 
features are naturally accounted for within our theory. However, we believe that the most 
remarkable and compelling evidence for our proposal comes from the observed coincidence 
of the burst activity with a large glitch. Moreover, the time evolution of the unabsorbed 
flux from AXP IE 2259+586 following the 2002 June outburst reported in Ref. jHO] can 
be explained naturally within our theory, but it is completely unreachable within the 
standard magnetar theory. So that we consider the June 18 SGR-like bursts from AXP 
IE 2259+586 the Rosetta Stone for our magnetar theory. 

The temporal decay of the flux during the burst activity displays a rapid initial decay 
which lasted about 1 days, followed by a more mild decay during the year following the 
onset of the burst activity. Indeed, the authors of Ref. [60J splitted the data into two 
segments, and fit each independently to a power law: 

F(t) ~ t ai , ai = - 4.8 ±0.5 , t < 1 days , 

F(t) ~ t a2 , a 2 = - 0.22 ± 0.01 , t > 1 days . ^ ' ' 

It is evident from Eq. ()5.18|) that the standard magnetar model is completely unable to 
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Figure 7: The time evolution of the unabsorbed flux from AXP IE 2259+586 following 
the 2002 June outburst. Data have been extracted from Fig. 13 in Ref. jHO]. Red dashed 
and dot-dashed lines are the phenomelogical power law fits, Eq. (|5.18|) . Blue continuous 
line is our light curve Eq. (J5.20|) . with parameters in Eq. (J5.23j) . 



reproduce the phenomelogical power law fit. On the other hand, even the phenomeno- 
logical parametrization cannot account for the time evolution of the flux. Indeed, if we 
assume the power law Eq. (j5.18|) for the decay of the flux, then we cannot explain why 
and how the source returns in its quiescent state with quiescent flux [60J: 

F Q ~ 1.53 lO- 11 ^^- . (5.19) 

cm z sec 

Note that adding the quiescent flux to the power law decay does not resolve the problem, 
for in that case the fits worst considerably. Our interpretation of the puzzling light 
curve displayed in Fig|7| is that AXP IE 2259+586 has undergone a giant burst at the 
glitch epoch, and soon after the pulsar has entered into a intense seismic burst activity. 
Accordingly, the flux can be written as: 

F{t) = F GB (t) + F sb if) + F Q , (5.20) 

where Fq is the quiescent flux, Eq. ()5.19|) . Fc B f) is the giant burst contribution to the flux 
given by Eq. ()5.6|) . and F SB it) is the seismic burst contribution given by Eq. (j5.17J) . Since 
during the first stage of the outburst there are no available data, we may parameterize 
the giant burst contribution as: 

VGB 

I t \ 1_ ''GB 

F GB {t) = F GB (0) 1 - — , < t < t GB , (5.21) 

V t GB J 
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while F$B(t) is given by: 



F SB {t) 



F S b(0) 

[l + mi)*™ 



i - 



ln(l + K X t) 
ln(l + K\t SB ) 



VSB 



< t < t 



SB 



(5.22) 



where toB and tgB are the dissipation time for giant and seismic bursts respectively. In 
Fig. we display our light curve Eq. (|5.2()jl with the following parameters: 



F GB (0) 
F SB (0) 



1.5 10" 



ergs 
cm 2 sec 



Vgb 



,, eras 
5.0 10~ n — ^— , V sb 

cm z sec 



0.828 , t GB ~ 0.91 days 
0.45 , t S B — 10 3 days , K\ 



0.20 days' 1 . 

(5.23) 



A few comments are in order. First, the agreement with data is rather good. Second, 
our efficiency exponent vgb is consistent with the values found in the giant bursts from 
SGR 1900+14 and SGR 0526-66, and cannot be justify in the standard magnetar model. 
On the other hand, quite consistently, we have tjsb < Vgb- Finally, we stress that from 
our interpretation of the light curve it follows that the onset of the intense seismic burst 
activity (KitsB ~ 200) did not allow a reliable estimation of which we predicted to be 
of order 10~ 2 . 

Interestingly enough, following the 2002 June outburst it was detected a infrared flux 
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Figure 8: The time evolution of the unabsorbed IR flux from AXP IE 2259+586 following 
the 2002 June outburst. Data have been extracted from Fig. 1 in Ref. |85j . Red dashed 
line is the phenomelogical power law fit t~°- 21±0 - 01 . Blue continuous line is our light curve 
Eq. (I5~m 
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changes correlated with the X-ray flux variability |85j . Since the observations began 
three days after the 2002 June outburst, according to our theory the infrared flux is 
parameterized as: 



t?IR(±\ F SB (0) 



V$ B 

r 1h(1 + ki*) " " 



;i + Kxt^ss 



ln(l + Kit SB ) 



1 — VSB 

+ Fq R , < t < t SB , (5.24) 



with the same parameters as in Eq. (|5.22j) . Indeed, assuming Fg B (0) ~ 9.5 10~ 15 crn -2 9 s s ec -i 
and Fq R ~ 3.3 10 15 cm -i 9 s s ec -i ; we found that our light curve Eq. ()5.24|) is in remarkable 
good agreement with data (see Fig. |HJ). The strong correlation between infrared and X- 
ray flux decay observed after the 2002 June outburst from AXP IE 2259+586 strongly 
suggest a physical link between the origin of both type of radiation. In particular, this 
unambiguously implies that the infrared flux originates from the magnetosphere. Indeed, 
in Appendix we argue that the soft emission from gamma-ray repeaters, anomalous X- 
ray pulsars, and isolated X-ray pulsars originates from thermal photons reprocessed by 
electrons trapped above the polar cups. On the other hand, we have already shown that 
the burst activity produces an increases of the thermal flux from the pulsar surface. So 
that, we see that in our theory the observed correlation between infrared and X-ray flux 
decays finds a natural explanation. 

5.2 SGR 1900+14 

Soon after the 1998, August 27 giant burst, the soft gamma repeater SGR 1900+14 entered 
a remarkable phase of activity. On August 29 an unusual burst from SGR 1900+14 
was detected jHS] which lasted for a long time ~ 10 3 sec. As discussed in Ref. [EE], on 
observational grounds it can be ruled out extended afterglow tails following ordinary 
bursts. Moreover, the standard magnetar model predicts faint transient afterglows on a 
time scale comparable to the duration of the bright X-ray emission of the burst peak. So 
that such mechanism cannot explain the long extended afterglow tail of August 29. 

In Figure M we display the flux decay after the August 29 burst. Data has been 
extracted from Ref. [SSI- m EE] the temporal behavior of the flux decay has been 
parameterized as a power law (red line in Fig. [SJ): 

F(t) = (89.16 ± 1.34) 10- 10 r(°-602±o.o25) . ( 5 . 25 ) 

cm 1 sec 

As already stressed, the phenomenological power law decay cannot explain the return of 
the source in its quiescent state where the flux is |87j : 

F Q = 0.96 ± 0.07 10~ n 6 \ gS . (5.26) 

cm 2 sec 

On the other hand, we may easily account for the observed flux decay by our light curve: 

ln(l + ^t) 



F(t) = m 
[) (1 + Kl t)> 



V 

1-1 



i , I . , + Fq , (5.27) 
ln(l + Kit di 

where Fq is fixed by Eq. (|5.26jl . Indeed, in Fig. El we compare our light curve Eq. ()5.27j) 
with observational data. The agreement is quite satisfying if we take: 

F(0) ~ 1.05 10" 9 6 \ gS , n ~ 0.5 , t dis ~ 1.2 10 3 sec , «i ~ 0.50 secT 1 . (5.28) 
cm 1 sec 



38 




time (sec) 



Figure 9: Flux evolution of the August 29 burst from SGR 1900+14- Data has been 
extracted from Fig. 4, panel (d), of Ref. [SHI- Red line is the phenomenological power law 
fit Eq. dHZHJ); blue line is our light curve Eqs. (fo~2"TJ) and fl£2EJ. 

The authors of Ref. |HH] have analyzed a large set of X-ray observations of SGR 
1900+14 in order to construct a more complete flux history. They found that the flux 
level was more than an order of magnitude brighter than the level during quiescence. This 
transient flux enhancement lasts about 40 days after the giant flare. Unlike the authors 
of Ref. [HE], which argued that this enhancement was an artifact of the August 27 flare, 
we believe that the flux history can be adequately described as seismic burst activity of 
the source. In Fig. El we report the flux light curve extracted from Fig. 2 of Ref. (HHJ 
together with their power law best fit. Again we find the the flux history is accounted for 
quite well by our light curve Eq. ()5.27j) with the following parameters: 

F(0) ~ 4.8 KT 8 6r / g , n ~ 0.55 , t dis ~ 200 days , k x ~ 2 10 3 days' 1 . (5.29) 
cm 2 sec 

The agreement between our light curve Eqs. (|5.27J) and ()5.29|) with the power law best fit 
is striking. Moreover we see that our curve deviates from the power law fit for t > 60 days 
tending to Fq at t = tdis- The authors of Ref. [HE] noted that extrapolating the fit to 
the August 27 X-ray light curve back toward the flare itself, one finds that the expected 
flux level lies below the ASM flux measurements (blue open points in Fig. 110)1 . Moreover, 
these authors observed that the discrepancy reduces somewhat when one pushes forward 
the reference epoch to about 14 minutes after the onset of the flare. However, unlike 
the authors of Ref. [88] we believe that the discrepancy is due to a true physical effect, 
namely the observed discrepancy from extrapolated light curve and ASM measurements 
is a clear indication that the surface luminosity increases after the burst activity. In 
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Figure 10: The time evolution of the unabsorbed flux from SGR 1900+14 following the 
1998 August outburst. Data has been extracted from Fig. 2 of Ref. [88J. Red line is the 
power law best fit F(t) ~ H°- 713±0 - 025 ) . Blu line is our light curve Eqs. (l5~2TI) and (fo~2"H|) . 



particular, soon after the August 27 giant flare we have seen in Sect. 14.31 that the surface 
temperature increases up to ~ 61 KeV and the surface luminosity reaches ~ 10 44 — . 

1 L J sec 

Almost all the deposited energy is dissipated within the dissipation time of the giant flare 
~ 400 sec. Nevertheless, it is natural to expect a more gradual afterglow where a small 
fraction of the energy deposited onto the star surface is gradually dissipated. As a matter 
of fact, we find that the observed level of luminosity Lx ~ 10 38 — (assuming a distance 
d = 10 kpc) at about 0.01 days since the August 27 giant flare, is consistent with the 
gradual afterglow scenario. 

On 2001 April 18 the soft gamma ray repeater SGR 1900+14 emitted an intermediate 
burst. The light curve of this event did not show any initial hard spike and was clearly 
spin-modulated. Moreover, the energetics appeared to be intermediate in the 40—700 KeV 
range, with a total emitted energy of about 1.9 10 42 ergs jHH]. In Fig. ITT1 we report the 
temporal behavior of the X-ray (2 — 10 KeV) flux from SGR 1900+14 in the aftermath 
of the 2001 April 18 flare. Data has been extracted from Fig. 2 of Ref. [EH!- The authors 
of Ref. |HH| attempted a simple power law function to the flux data: 

F(t) ~ t~ a + K , (5.30) 

where the constant K should take care of the quiescent luminosity. Indeed, the authors 
of Ref. H3 fitting Eq. ijOty to the data found: 

a = 0.89(6) , K = 0.78(5) 10~ n 6T f S . (5.31) 

cm 1 sec 
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Figure 11: Temporal behavior of the X-ray flux from SGR 1900+14 in the aftermath of 
the 2001 April 18 flare. Data have been extracted from Fig. 2 of Ref. [89]. Red line is 
the power law best fit Eqs. (|5.30|) and (|5.31jl . Blu line is our light curve Eqs. ()5.27j) and 

As it is evident from Fig. ^2 the power law globally fits the data quite nicely. However, 
the reduced x 2 turns out to be in excess to 3, mainly due to the bump occurring in the 
light curve at t ~ 10 5 sec jSHJ- Indeed, after excluding the bump they get a good fit with 
X 2 /dof ~ 1 without an appreciable variation of the fit parameters [EH]- However, there 
is still a problem with the phenomelogical power law decay of the flux. As a matter of 
fact, Eq. (J5.31)) shows that the power law fit underestimates the quiescent luminosity. In 
our opinion this confirms that the phenomenological power law decay of the flux is not 
adeguate to describe the time variation of the flux. On the other hand, we find that our 
light curve Eq. (|5.27|) . with quiescent luminosity fixed to the observed value Eq. ()5.26|h 
furnishes a rather good description of the flux decay once the parameters are given by: 

F(0) ~ 2.6 10~ 7 6r f S , r) ~ 0.68 , t dis ~ 3 10 6 sec , k x ~ 0.25 sec~ l . (5.32) 
cm 1 sec 

Within our interpretation, Eq. (|5.32|) shows that the flux decay in the aftermath of the 
April 18 flare is characterized by a very large seismic burst activity (nitdis ~ 10 6 ), which 
lasts for about 10 6 sec. So that the bump in the flux at t ~ 10 5 sec is naturally explained 
as fluctuations in the intensity of the seismic bursts. 

In Ref. IHIIj it is reported the spectral evolution and temporal decay of the X-ray tail of 
a burst from SGR 1900+14 recorded on 2001 April 28, 10 days after the intense April 18 
event. In Fig.^Jwe display the temporal decay of the flux from the 2001 April 28 burst in 
the energy band 2 — 20 KeV. The data has been extracted from Fig. 5 of Ref. [HHl ■ These 
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Figure 12: The temporal decay of the flux from the 2001 April 28 burst from SGR 
1900+14 in the energy band 2 — 20 KeV. The data has been extracted from Fig. 5 of 
Ref. PHI- Red continuous line is the phenomenological best-fit power law times exponential 
function adopted in Ref. jUHl to describe data, Eqs. ()5.33|) and (|5.34j) . Blue continuous 
line is our light curve Eqs. (|5.27j) and ([5.35)1 . 



authors attempted several functional forms to fit the decay of the flux. They reported that 
the decay was equally well fitted by a either power law times exponential or broken power 
law. We stress that both fits are phenomenological parametrization of the observational 
data, and that both fits are unable to recover the quiescent flux. For definitiveness, we 
shall compare our light curve with the power law times exponential fit: 

F(t) ~ t~ a exp(--) , (5.33) 

T 

The best fit to the temporal decay of the flux from the 2001 April 28 burst in the energy 
band 2-20 KeV gives |H]: 

a = 0.68 ± 0.04 , r = 5 ± 1 10 3 sec . (5.34) 

In FigureEDwe compare the phenomenological best fit Eqs. (j5.33J) and (|5.34j) with our light 
curve Eq. (|5.27|) , where the quiescent luminosity is fixed to the observed value Eq. (|5.26|) , 
and the parameters are given by: 

F(0) ~ 1.4 10" 9 6r 9S , T) ~ 0.5 , t dis ~ 5.5 10 3 sec , k x ~ 0.06 sec' 1 . (5.35) 
cm 2 sec 

Again, we see that our light curve gives a quite satisfying description of the flux decay. 
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Figure 13: The temporal decay of the flux from the 1998 August 29 burst from SGR 
1900+14 in the energy band 2 — 20 KeV. The data has been extracted from Fig. 10 of 
Ref. 90j. Red continuous line is the phenomenological best-fit power law times exponential 
function adopted in Ref. [HOI to describe data, Eqs. (|5.33|) and (|5.36|) . Blue continuous 
line is our light curve Eqs. (|5.27j) and (|5.37|) . 

Interestingly, the authors of Ref.[9*U] analyzed with the same techniques the 1998 
August 29 burst from SGR 1900+14- The fit to the decay of the flux in the energy band 
2-20 KeV with Eq. (I5~531) resulted in: 

a = 0.510 ± 0.008 , r = 2.9 ± 0.2 10 3 sec . (5.36) 

Even in this case our light curve is able to follow the time decay of the flux in a satisfying 
way. In Fig. EH we compare our light curve with parameters: 

F(0) ~ 5.0 10~ 9 6r / g , V - 0.45 , t dis ~ 7.0 10 3 sec , m ~ 0.25 sec' 1 , (5.37) 
cm 2 sec 

and the phenomelogical power law times exponential fit Eqs. (j5.33|) and ()5.36|) . 

Let us consider, finally, the light curve for the intermediate burst from SGR 1900+14 
occurred on 2001 July 2 [§!]. In Figure IT41 we display the time decay of the flux after 
the July 2 burst. The data have been extracted from Fig. 1 of Ref. by binning the 
light curve histogram. The displayed errors are our estimate, so that the data should be 
considered as purely indicative of the decay of the flux. We find that Fig. 1 of Ref. [01] is 
very suggestive, for it seem to indicate that the burst results from several small bursts, 
i.e. according to our theory the burst is a seismic burst. As a consequence we try the fit 
with our light curve Eq. ([5.27)1 . In this case the quiescent flux has been fixed to Fq ~ 0, 
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Figure 14: Time history of the 2001 July 2 burst from SGR 1900+14 in the energy band 
7 — 100 KeV as observed by FREGATE. The data has been extracted from Fig. 1 of 
Ref. [HI] by binning the light curve histogram. Blue continuous curve is our light curve 
Eqs. (I5~?7f) and (ET3SI) . 



for the observational data has been taken in the energy range 7 — 100 KeV where the 
quiescent flux is very small. Attempting the fit to the data we find: 

coiiTvbs 

F(0) ~ 4.07 10 4 , T] ~ 0.36 , t dis ~ 40 sec , «i ~ 5.0 sec -1 . (5.38) 

sec 

The resulting light curve is displayed in Fig. HU The peculiarity of this burst resides in 
the fact that the burst activity terminates suddenly at t ~ 4 sec well before the natural 
end at t^ s ~ 40 sec. 



5.3 SGR 1627-41 

SGR 1627-41 was discovered with the Burst And Transient Source Experiment (BATSE) 
on the Compton Gamma- Ray Observatory (CGRO) in June 1998 [2211231 when it emitted 
over 100 bursts within an interval of 6 weeks. The authors of Ref |2j presented the results 
of the monitoring of the flux decay of the X-ray counterpart of SGR 1627-41 spanning an 
interval of roughly five years. Moreover, these authors attempted to understand the three 
year monotonic decline of SGR 1627-41 as cooling after a single deep crustal heating 
event coinciding with the burst activity in 1998 within the standard magnetar model. 
They assumed an initial energy injection to the crust of the order of 10 44 ergs. However, 
it must be pointed out that this assumption is highly unrealistic, for the estimate of 
the total energy released in bursts during the activation of SGR 1627-41 range between 
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Figure 15: Time decay of the flux from SGR 1627-41- The data has been taken from 
Table 1 of Ref. j^lj. Blue continuous line is our best fitted light curve Eqs. ([5.27]) . ()5.39|) 
and (EH- 



4 10 — 2 10 43 ergs. In addition, since gamma rays was not detected, they assume that 
the conversion efficiency of the total energy released during the activation into soft gamma 
rays were considerable less than 100 %. They also assumed that the core temperature is 
low, i.e. the core cools via the direct URCA process [30]. Notwithstanding these rather 
ad hoc assumptions, the authors of Ref. [53] was unable to explain the March 2003 data 
point, which clearly showed that the flux did not decay further (see Fig. I15|) . In other 
words, the levelling of the flux during the third year followed by its sharp decline is a 
feature that is challenging that standard magnetar model based on neutron stars, and 
that beg for an explanation within that model. On the other hand, we now show that 
the peculiar SGR 1627-41 light curve find a natural interpretation within our theory. In 
Fig. E3 we display the time decay of the flux. The data has been taken from Table 1 
of Ref. [HI]. In this case we are able to best fit our light curve Eq. (J5.27|) to available 
data. Since the number of observations is rather low, to get a sensible fit we have fixed 
the dissipation time to 1200 days and the quiescent luminosity to the levelling value at 
t > 1200 days: 

F Q ~ 2.7 10- 13 e [ 9S , t dis ~ 1200 days . (5.39) 
cm* sec 

The best fit of our light curve to data gives: 

F(0) = 0.83(11) 10~ n 6r / g , r] = 0.25(8) , m = 0.04(1) days' 1 . (5.40) 

cm* sec 

with a reduced \ 2 ~ 1- From Fig. where we compare our best fitted light curve 
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with data, we see that our theory allow a quite satisfying description of the three year 
monotonic decline of the flux. 

5.4 SGR 1806-20 
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Figure 16: Time history of the second component after the 2004 December 27 giant burst 
from SGR 1806-20. The data has been extracted from Fig. 5 of Ref. jHS] by binning the 
light curve histogram. Blue continuous curve is our light curve Eqs. ([5.27)1 . (|5.42j) and 

(JOB)- 

SGR 1806-20 entered an active phase in 2003, culminating in a gigantic flare on 2004 
December 27, with energy greatly exceeding that of all previous events. In Figure 3 of 
Ref. [80 it is reported the time history of the flux averaged over the rotation period of the 
pulsar soon after the giant flare. These authors fitted the light curve within the standard 
magnetar model based on the evaporation of a fireball formed after the giant burst and 
trapped onto the stellar surface, Eq. ([5.7)1 . The fit of the rotation smoothed curve to the 
fireball function gives 80J: 

a = 0.606 ± 0.003 , t evap = 382 ± 3 sec . (5.41) 

However, from Figure 3 of Ref. jHOl it is evident that fireball light curve underestimate 
the luminosity for t < 30 sec. Thus we see that the light curve can be better accounted 
for by our light curve Eq. ([5.6)1 with tb rea k — 30 sec, quite close to the values found for 
the giant bursts from SGR 1900+14 and SGR 0526-66. A more precise determination 
of the parameters of our light curve, however, must await for more precise data in the 
initial phase of the afterglow. Instead, in the present Section we discuss the light curve of 
a second, separate component after the giant burst reported in Ref. The authors of 
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Ref. [95] found evidence for a separate component in the light curve starting at t ~ 400 sec 
from the onset of the giant burst, forming a peak at t ~ 600 sec and ending at t ~ 3000 sec 
(see Fig. 5 of Ref. [EH])- As already discussed, in our theory it is expected that there is 
a intense seismic burst activity following a giant burst. In Figure IT^l we display the flux 
history starting from the giant flare. We show a few point of the second component 
extracted from Fig. 5 of Ref. [HHj by binning the light curve histogram. The displayed 
errors are our estimate, so that the data should be considered as purely indicative of the 
decay of the flux. We fit the data with our light curve Eqs. f)5.27|) assuming: 

, count 

Fq ~ 8.796 10 , t dis = t end - t star t , t sta rt — 625 SCC , t end ~ 3000 SCC . 

sec 

(5.42) 

The best fit of our light curve to data gives: 

F(0) = 0.074 10 4 , r) = 0.18 , m = 0.10 sec' 1 . (5.43) 

sec 

Indeed, in Fig. El we compare our light curve with data and find that our theory allow a 
quite satisfying description of the time history of the flux. 



6 CONCLUSIONS 



Let us summarize the main results of the present paper. We have discussed p-stars en- 
dowed with super strong dipolar magnetic field. We found a well defined criterion to 
distinguish rotation powered pulsars from magnetic powered pulsars (magnetars). We 
showed that glitches, that in our magnetars are triggered by magnetic dissipative effects 
in the inner core, explain both the quiescent emission and bursts from soft gamma-ray 
repeaters and anomalous X-ray pulsars. In particular, we were able to account for the 
braking glitch from SGR 1900+14 an d the normal glitch from AXP IE 2259+586 fol- 
lowing a giant burst. We accounted for the observed puzzling anti correlation between 
hardness ratio and intensity. Within our magnetar theory we were able to account quan- 
titatively for light curves for both gamma-ray repeaters and anomalous X-ray pulsars. 
In particular we explained the light curve after the June 18, 2002 giant burst from AXP 
IE 2259+586. Finally, in Appendix we discussed the origin of the soft emission from soft 
gamma-ray repeaters, anomalous X-ray pulsars, isolated X-ray pulsars. 
We believe that anomalous X-ray pulsars and soft gamma-ray repeaters are two class of 
intriguing objects that are challenging the standard paradigm based on neutron stars. 
In the present paper we convincingly argued that the standard magnetar theory is com- 
pletely unable to account for the observational properties of anomalous X-ray pulsars and 
soft gamma-ray repeaters. On the other hand, we feel that the ability of our p-star theory 
to reach a complete understanding of several observational features of soft gamma-ray 
repeaters and anomalous X-ray pulsars strongly supports our proposal for a drastic revi- 
sion of the standard paradigm of relativistic astrophysics. 

Let us conclude by briefly addressing the theoretical foundation of our theory. As a 
matter of fact, our proposal for p-stars stems from recent numerical lattice results in 
QCD jHElinZI, which suggested that the gauge system gets deconfined in strong enough 
chromomagnetic field. This leads us to consider the new class of compact quark stars 
made of almost massless deconfined up and down quarks immersed in a chromomagnetic 
field in /^-equilibrium. Our previous studies showed that these compact stars are more 
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stable than neutron stars whatever the value of the chromomagnetic condensate. This 
remarkable result indicates that the true ground state of QCD in strong enough gravita- 
tional field is not realized by hadronic matter, but by p-matter. In other words, the final 
collapse of an evolved massive star leads inevitably to the formation of a p-star. 

A ORIGIN OF THE SOFT EMISSION IN X-RAY PULSARS 

A number of anomalous X-ray pulsars have recently been detected in the optical-infrared 
wavelengths (for a recent review, see Ref. |HE|)- Since pulsar surface emission cannot 
account for the observed soft emission spectra, the emission must be magnetospheric 
in origin. This is clearly demonstrated by the correlation of the infrared flux with the 
bursting activity recently observed for the anomalous X-ray pulsar AXP IE 2259+586 [85J 
(see Sect. 15.1(1 . Remarkably, there is compelling evidence for an excess of the flux in 
the infrared band in anomalous X-ray pulsars with respect to the thermal component 
extrapolated from X-ray data. Up to now the standard magnetar model, is unable to 
account for the observed infrared emission or variability. In any case, irrespective to the 
actual mechanism responsible for the soft emission there is no doubt that the optical- 
infrared emission originates in the magnetosphere. Looking at the broad band energy 
spectrum of anomalous X-ray pulsars one realizes quickly that the spectra are amazingly 
similar to the ones of isolated X-ray pulsars, like RXJ 1856.5-3754 or RXJ 0720.4-3125. 
This strongly suggest that it must exist some natural mechanism capable to generate the 
soft tail of the spectrum of isolated pulsar, anomalous X-ray pulsars and soft gamma ray 
repeaters. In Refs. EI] we already advanced the proposal that the faint emissions from 
RXJ 1856.5-3754 an d RXJ 0720.4-3125 originate in the magnetosphere from synchrotron 
radiation emitted by electrons. That proposal was based on the observational fact that 
the faint emission can be parameterized quite well by a non thermal power law. However, 
we did not address the problem of the physical mechanism which is at the heart of the 
electron energy spectrum needed to generate the power law emission. In this Appendix 
we shall discuss a fair natural mechanism which is able to explain the faint emission for 
isolated pulsars as well as anomalous X-ray pulsars and soft gamma ray repeaters. In 
our mechanism the power law emission in the infrared-optical band is due to thermal 
radiation reprocessed in the magnetosphere by electrons trapped near the magnetic polar 
cups. 

To start with, let us consider the motion of a charged particle in the pulsar magnetic 
dipolar field. We are assuming the presence of neutral plasma formed by electrons and 
protons with number densities n e = n p . Now, as is well known, for strong enough magnetic 
fields the plasma near the stellar surface will be channelled toward the magnetic pole. So 
that we are led to consider the motion of charged particles which are drifting toward the 
magnetic polar cups. Further, we may consider a polar region with area small enough 
such that the dipolar magnetic field depends only on the distance from the surface. If the 
magnetic polar axis is in direction z, then charged particles moving towards the stellar 
surface will feel an almost uniform magnetic field having z-component B(z). Thus we are 
led to consider the motion of charge particle in the magnetic field: 

B{z) = -B s ^j , z > R , (A.l) 
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where, for the sake of definitiveness we are considering the north magnetic pole where 
the magnetic field is entering into the stellar surface. Let us consider, firstly, electron 
with charge — e and mass m e . The electron wave function ip(x,y,z) can be obtained by 
solving the Schrodinger equation in presence of the magnetic field B(z). We are interested 
in the physical problem where an electron, starting from a distance zq from the star, is 
approaching the stellar surface. Usually, it is assumed that the magnetic field is almost 
constant, so that our problem reduces to the well known motion in an uniform magnetic 
field. In this case one gets: 



^{x,y,z) = exp(-ip z z) <f) n (x,y) 



with energy eigenvalues: 



PI 



■n,p z 



2m f 



eB s 1,1, 
n + - ± - 

m e K 2 2' 



n 



0,1, 



(A.2) 



(A.3) 



However, if z ^> R the uniform field approximation is no longer valid, but the weakly 
varying field is more appropriate. In this case we may write: 

ip(x,y,z) = ((z) <fi n (x,y) , (A.4) 

where <j) n (x,y) is the solution of the Schrodinger equation in the weakly varying B(z): 



2m e 



|^ + e 2 5W + 2iB(z)^- 
ox z ox 



dy 2 
1, 



<j> n (x,y) = e n (z)(j) n (x,y) 



eB{z) . 1 
e n (z) = — ^ n + - ± - 
m P 2 2 



n 



0,1, 



So that, writing the total energy of electrons drifting from zq toward the star as: 

E = e n (z Q ) + e drift , 
the wave function ((z) satisfies the effective Schrodinger equation: 

1 



2m„ 



where: 



VJz) 



eBs R 3 



m P Zn 



■J? + Vs{z) 



n H — ± - 

2 2' 



C(z) = e drift C(z) 



2(1 



R < Z < Zq 



(A.5) 
(A.6) 

(A.7) 
(A.8) 

(A.9) 



Inspection of Eq. (JA.9j) shows that, as long as the total momentum is not exactly parallel 
to the magnetic field, electrons will feel a repulsive barrier which at the stellar surface is: 



eB< 



11.6 KeV B 12 , B 



Be 



12 



(A.10) 



m e 10 12 Gauss ' 

at least. This effect is the quantum mechanical counterpart of the well known fact that 
classical charges moving towards regions with increasing magnetic fields are subject to a 
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Figure 17: Blue line is the effective potential Eq. (|A.13|) for electrons drifting toward the 
stellar surface in the polar cup regions. The red line is the harmonic approximation to 
the potential. 



repulsive force. On the other hand, when we consider protons the minimal height of the 
magnetic barrier is: 

V = — ~6.3B 12 eV . (A.ll) 



m r 



In other words, protons drifting towards the surface almost do not feel any barrier. As a 
consequence, there is an accumulation of positive charge on the surface. Note that on the 
surface of the star there is a positively charged layer which is able to support a thin crust 
of ordinary matter. Thus, protons do not dissolve into the star, but they are trapped 
in the atmosphere of electrons which extends over a distance 5 ~ 10 3 fermis beyond the 
edge. Let n be the number density of trapped protons, then on the surface of the star 
there is a surface charge density a ~ 5 n which, in turn, gives rise to an uniform electric 
field. It follows that electrons that are moving toward the star are subject to both the 
repulsive magnetic potential and the attractive electric potential. So that the effective 
potential is: 



V{z) 



eB s R 3 



1 



1. 



(n + - ± - 



1 



+ 4ne 2 n5z, R < z < z . (A.12) 



In the case of minimal magnetic barrier, using Eq. (jA.lO)) we rewrite Eq. (jA.12j) as: 



R 3 

V{z) ~ 11.6 KeV B 12 —» 



- _ i 



+ 0.23 KeV m 3 - , 1 
K 



Z Zn 

< — < — 
- R - R 



(A.13) 
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where n% 3 = 10l3 " m _ 3 . In Fig. [T7| we display the effective potential Eq. (|A.13|) assuming 
z ~ 10 R. We see that the effective potential V(z) displays a minimum at z = z where 
the repulsive magnetic force is balanced by the attractive electric force. Then, electron are 
trapped above the polar cup at a distance of order z. To determine the energy spectrum 
we need to solve the Schrodinger equation with the effective potential: 



1 



1m „ 



d 



2 



+ V(z) 



C(z) = e dri f t ((z) . (A.14) 



We may adopt the harmonic approximation to the potential by expanding around z. A 
straightforward calculation gives: 



1 13 



z ~ 3.5 R B* 2 n 13 J , V = V(z) ~ 1.08 KeV Bf 2 nf 3 . (A.15) 

Moreover: | 

R 2 V"(z) ~ 0.261 KeV B~} nf 3 . (A. 16) 

Note that, as long as zq 3> R, z, V and V (z) do not depend on z . In this approximation 
((z) satisfies the harmonic oscillator equation centered at z with frequency: 

uj m ~ 0.46 10~ 12 eV nf 3 . (A.17) 
Thus, we find for the drift energy the quasi continuum spectrum: 

£drift,j = V + uo m (j + -) , j > . (A.18) 

Let us pause to briefly summarize our results. Our quantum mechanical treatment of 
charges which are drifting toward the star in the weakly varying dipolar magnetic field 
has shown that electrons feel a huge magnetic barrier. On the other hand, the magnetic 
barrier is reduced by a factor ^ for protons. As a consequence, protons are free to reach 
the stellar surface, where they are trapped in the electron atmosphere, while electrons are 
repelled into the magnetosphere. The resulting charge separation produces an electric field 
which trapped electron at a distance ~ z. Obviously, the neutrality of the plasma implies 
that the number densities of trapped electrons and protons are equal. The resulting 
picture is very interesting. We have electrons which oscillate around z along the magnetic 
field. On the other hand, these electrons feel the magnetic field B(z), which now is truly 
almost constant. So that, if u B is the cyclotron frequency at z, then using Eq. (|A.15|) we 
get: 

eB(z) i 2 

u B = — — ~ 0.268 KeV Bf 2 nf 3 . (A. 19) 

Then, the energy spectrum of these electrons is: 

E nd = V + oj b n + ou m (j + ~) , n , j > . (A.20) 

In Figure^] we illustrate schematically the spectrum Eq. (|A.20J) . According to Eq. (|A.20J) 

the spectrum comprises discrete Landau levels with cyclotron frequency uj b and an almost 
continuum associated to the drifting motion. As discussed below, electrons perform ra- 
diative transition between Landau levels. These transitions are responsible for absorption 
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Figure 18: Energy spectrum for electrons trapped in the magnet osphere near the polar 
cups. Heavy lines are the Landau levels with cyclotron frequency ujb- Light lines are the 
quasi continuum spectrum Eq. (jA.18|) . Red and blue arrows indicate possible radiative 
transitions. 

and spectral features observed in isolated X-ray pulsars, anomalous X-ray pulsars and 
soft gamma ray repeaters. Moreover, we see that there are also transitions where elec- 
trons absorb thermal photons with frequency uj ~ ojb and emit photons with frequencies 
uj <C uj. These radiative transitions give rise to the soft spectrum. Before addressing the 
problem of the soft spectrum, let us discuss the puzzling absorption features detected in 
the isolated X-ray pulsar IE 1207.4-5209 [H]. The spectrum of IE 1207.4-5209 shows 
three distinct features, regularly spaced at 0.7 KeV, 1.4 KeV and 2.1 KeV, plus possibly 
a fourth at 2.1 KeV. These features vary in phase with the star rotation. Indeed, it 
turns out that the X-ray source pulsation is largely due to the phase variation of the 
lines with the pulsar rotation. The most natural and logical explanation for the observed 
features is cyclotron resonant absorption. In this case the fundamental cyclotron fre- 
quency is 0.7 KeV. Thus, the inferred magnetic field is ~ 8 10 10 Gauss for electrons and 
~ 2 10 14 Gauss for protons. However, the magnetic field inferred from the spin parame- 
ters turns out to Bs — 2.4 10 12 Gauss. Obviously, as noted in Ref. [§H| electron cyclotron 
scattering at a distance 3 — 4 R above the pulsar surface would fit all the observations. 
Remarkably, using the inferred magnetic field Bs and Eq. (jA.19|) we find: 

cob ^ 0.7 KeV , B s ,i2 ^ 2.4 n 13 ~ 2.69 . (A.21) 

Moreover from Eq. (|A.15|) we get z ~ 3.4 R, in perfect accord with observations. 
Finally, we stress that the presence of the almost continuum spectrum associated to 
the drifting motion along the magnetic field leads naturally to rather broader spectral 
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features in accord with several observations. Even though it is beyond the aim of the 
present Section an accurate comparison with available data, we would like to stress that 
our proposal is in gratifying qualitative agreement with observations. Indeed, absorption 
features have been found in the thermal emission of several isolated X-ray pulsars which 
range in 0.2 — 0.7 KeV. The thermal spectrum of isolated X-ray pulsar is a blackbody 
with typical temperature T ~ O.lKeV. Now Eq. (jA.19|) shows that for typical pulsar 
magnetic field B12 ~ 1 the electrons trapped at ~ z are able to absorb thermal photons 
in the observed range. Note that for magnetars the typical blackbody temperature is 
T ~ 1.0 KeV (the relevant temperature is the polar cup blackbody temperature), and 
B12 ~ 10 2 . So that Eq. (jA.19|) gives u B ~ 0.85 KeV, indicating that also for magnetars 
the trapped electrons may efficiently absorb the thermal photons from the polar cups. 
Let us discuss, now, the radiative transitions from trapped electrons. To evaluate the rate 
of transitions we need to evaluate: 



< n ,j I exp(— ik ■ r) I ■ V |n, j > 



where: 



ki > = CM) <t>n{x,y) 



I ■ k 







(A.22) 



(A.23) 



k and / being the photon momentum and polarization respectively. Without loss in 
generality, we may write: 



k = (k x ,0,h) , I = (0,1,0) , 



(A.24) 



so that: 



where: 



— < n ,j I exp(— ik ■ r)l ■ V In, j > 



Obviously we have also: 
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(A.26) 
(A.27) 
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E n',j' ~ E n,j 



u B (n - n) + uj m (j - j) 



(A.28) 



We are interested in processes where trapped electrons absorb thermal photons with 
energy ~ uj b (red lines in FiglTSj) and emits photons with \k\ -C ujb (blue lines in Fig.fTHj). 
Since u m ~ 10 -12 eV, from Eq. (jA.28|) it follows that Aj ^> 1. Then, the transitions 
j — > j cannot be induced by the electromagnetic field. Indeed, these transitions are 
induced by thermal collisions. To see this, we note that for j ^> 1 the wave function ((z) 
is quasi classical: 

CO) ~ exp[-ip(z)z] , (A.29) 

where p(z) is the quasi classical momentum. Now, Eq. (|A.27J) means that J-i ■ is different 
from zero if the momentum is conserved. However, for vastly different quasi momenta 
the overlap of the quasi classical wave function is very small. So that we are lead to the 
conclusion that | 1 <C \k±\. As a consequence we see that the transitions must be induced 
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by collisions, so that we may safely assume \Jj'j\ — 1. To evaluate I n r n we need the wave 
functions (f> n (x, y): 



exp(-tp x x) 



r exp[ 
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m e uj B 
(A.30) 

l,n — 0. A standard calculations gives: 
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so that: 
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Note that in our approximation \k\ — k \k±\ <C 0Jb- Thus, we have: 



\A 



if I 



ie 



< n ,j | exp(— ik ■ r) I -V \n,j > 



e 2 uj B 
2m„ 



2m e u B 

(A.32) 



(A.33) 



It is clear that Eq. (jA.33|) leads to a power law emission flux with the cutoff k uj b . 
Moreover, to get the flux we need to assume some energy distribution for the trapped 
electrons, which, however, goes beyond the aim of the present paper. Nevertheless, there 
is a general aspect which is worthwhile to stress. We have seen that the soft spectrum in 
isolated pulsar can be understood as thermal photons reprocessed by electrons trapped in 
the magnetosphere above the polar cups. These electrons absorb thermal photons with 
frequency ~ u B and emit photons with frequencies <^ uj b . So that the number of these 
electrons is proportional to the number of thermal photons with energy ~ uo B . We know 
that uo B ~ T, where T is the blackbody temperature. As a consequence, recalling that 
the number density of thermal photons scale as ~ T 3 , we have that the soft fluxes in 
magnetars should be about a factor {T magne tar/T pu i s ar)^ ~ 10 3 greater than the one in 
isolated X-ray pulsars, in fair agreement with observations. 

Let us conclude this Appendix by roughly estimating the soft flux F^. The probability 
for emission of a photon in the frequency range u, uj + du is: 



dP 

dt duj 
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5[uj - {E x , - Eoj)] 



(A.34) 



where we used ou — \k\ — k ^ \k±\, and the delta function ensures the energy conservation 
Eq. (|A.28|) . Summing over the degenerate final states of the almost continuum spectrum, 
and using du = uj m dj , we get: 
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1 uj e uj b 
2ix uj m 2m P 



(A.35) 



Finally, the spectral flux is obtained multiplying by the photon energy uj and by the 
number of active electrons. Let n act be the number density of active electrons, we have: 



1 uj 2 e 2 u B 
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2n uj m 2m e 
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(A.36) 
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where V is the volume of the emitting region. Note that F u has a Rayleigh- Jeans power 
law form with an upper cutoff which, however, cannot be easily estimated without a 
precise knowledge of electron energy distribution. Using Eqs. (|A.17|) and (jA.19|) . we 
rewrite Eq. (jA.36|) as: 

F u ~ 8.2 10- 6 (-^-) 2 Bf 2 ni n act V . (A.37) 

sec Hz lev 

The number density of active electrons is not easily estimated without the knowledge 
of the energy distribution of electrons trapped into the magnet osphere. In general, n act 
depends on the number of thermal photons with energy u ~ ujb- Nevertheless, we may 
estimate the needed number density n act by comparing with observed soft fluxes. For 
typical isolated X-ray pulsar with Z? 12 ~ 1 and n 13 ~ 1 we infer: 

F u ~ 10 13 6T9 T S T , u ~ leV . (A.38) 

sec Hz 

So that, assuming for the emitting volume the reasonable value V ~ 10 10 cm 3 , from 
Eq. (TOTI) we get: 

n act ~ 10 8 cm" 3 ~ I0~ 5 n . (A.39) 

It is interesting to compare Eq. (|A.39)1 with the number density of active electrons obtained 
assuming that trapped electrons have an uniform distribution. In this case, observing that 
the probability for transitions from the n = to n = 1 Landau levels is given by Eq. (lA~32l) 
with k± ~ u B , we get: 

"--^7"- < A - 40 > 

So that, for typical isolated X-ray pulsars we have: 

n act ~ 2.4 10~ 5 n , (A.41) 
that, indeed, is in reasonable agreement with our estimate Eq. (jA.39|) . 
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